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Abstract 

, We present in detail the calculation of the 0{oP s ) inclusive total cross section for the process 

> : 

t^- ■ PP ~ ¥ ^th, * n ^ ne Standard Model, at the CERN Large Hadron Collider with center-of-mass energy 

OO ' _ 

y / s^"=14 TeV. The calculation is based on the complete set of virtual and real 0(a s ) corrections 

m : 

to the parton level processes qq — * tth and gg — > tth, as well as the tree level processes (q, q)g — ► 

m ■ 

tth + (q,q). The virtual corrections involve the computation of pentagon diagrams with several 
internal and external massive particles, first encountered in this process. The real corrections 



Oh 



are computed using both the single and the two cutoff phase space slicing method. The next- 
to-leading order QCD corrections significantly reduce the renormalization and factorization scale 
dependence of the Born cross section and moderately increase the Born cross section for values of 
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the renormalization and factorization scales above mt- 
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I. INTRODUCTION 



One of the critical goals of present and future colliders is the study of the electroweak 
symmetry breaking mechanism and the origin of fermion masses. If the introduction of one 
or more Higgs fields is responsible for the breaking of the electroweak symmetry and for the 
generation of fermion masses, then one Higgs boson should be relatively light. The present 
lower bounds on the Higgs boson mass from direct searches at LEP2 are > 114.4 GeV 
(at 95% CL) [l| for the Standard Model (SM) Higgs boson (h), and M h o > 91.0 GeV and 
M A o > 91.9 GeV (at 95% CL, 0.5 < tan/3 < 2.4 excluded) Q for the light scalar (h°) and 
pseudoscalar (A ) Higgs bosons of the minimal supersymmetric standard model (MSSM). 
At the same time, global SM fits to electroweak precision data imply < 211 GeV (at 95% 
CL) P], while the MSSM requires the existence of a scalar Higgs boson lighter than about 
130 GeV. The possibility of a Higgs boson discovery in the mass range near 115-130 GeV 
thus seems increasingly likely. 

The associated production of a Higgs boson with a tt pair can play a very important role at 
hadron colliders as has been suggested by many studies over the past several years jj, ^, 0, E | . 
In particular, it is an important discovery channel for a SM-like Higgs boson at the LHC if 
Mh < 130 GeV ^,y,y,ll3- Although the event rate is small, the signature is quite distinctive. 
Given the statistics expected at the LHC, pp — ► tih, with h — > bb,T + r~ ,W + W~ ,77 will 
also be instrumental to the determination of the couplings of a discovered Higgs boson, and 
will in particul ar g ive the only handle on a direct measurement of the top quark Yukawa 

coupling 0, H HQ, Q . 

The total cross section for pp — > tih has been known at tree-level, i.e. at leading order 
(LO) of QCD, for quite some time 15, Next-to-leading order (NLO) QCD corrections 



are crucial in order to reduce the dependence of the cross section on the renormalization 
and factorization scales. The calculation of the total cross section for pp — > tih to 0(a^) 



has been performed by the Authors of Refs. 



17L Il8l | and by our group. The results of the 



two independent calculations have been compared and they are in very good agreement. In 
Ref. [l9j |. we presented our first numerical results for the total inclusive NLO QCD cross 
section for pp — > tih at the LHC center of mass energy, ■ s fs~^= 14 TeV. Here we provide a 
detailed description of the calculation. 

At the LHC center-of-mass energy, the dominant subprocess for tih production is gg — > 
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tih, but the other subprocesses, qq — > tih and (q,q)g — > it/i + (?,?), which contribute to 
the cross section at C?(a^), cannot be neglected and are included in this calculation. The 
NLO QCD corrections to the qq — > tih subprocess constitute a gauge i nvariant subset of the 
entire NLO QCD calculation and have been presented in Refs. 2C| |2l| to which we refer for 
a thorough discussion of the results. Here we concentrate on a detailed description of the 
calculation of the 0(a s ) corrections to the gg — > tih subprocess. The Feynman diagrams 
contributing to gg — > tih at lowest order are shown in Fig. ^ while the 0(a s ) virtual and 
real corrections are given in Figs. |2E and Fig. ® respectively. 

The main challenge in the calculation of the 0(a s ) virtual corrections comes from the 
presence of pentagon diagrams with several massive external and internal particles. The 
pentagon scalar and tensor Feynman integrals originating from these diagrams present either 
analytical (scalar) or numerical (tensor) challenges. We have calculated the pentagon scalar 
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integrals as linear combinations of scalar box integrals using the method of Ref. 
and cross checked them using the techniques of Ref. (2^]. Pentagon tensor integrals have 
been calculated and cross checked in two ways: numerically, by isolating the numerical 
instabilities and extrapolating from the numerically safe to the numerically unsafe region 
using various methods; and analytically, by reducing them to a numerically stable form. 
The real corrections have been computed using the phase space slicing method, in both 
the double (for a review see, e.g. |25j ) and sing 

le 0HQ 

cutoff approaches. Together 

with the corresponding qq — ► tih calculation [20, |2l[ , this is the first application of the single 



cutoff phase space slicing method to a cross section involving more than one massive particle 
in the final state and agreement between the two cutoff and the single cutoff approaches is 
a strong check of the calculation. 

The outline of our paper is as follows. In Section |H] we summarize the general structure 
of the NLO cross section for pp — > tih. In Section IIHI we briefly review the case of the 
LO cross section for pp — > tih, introducing some fundamental notation. We proceed in 
Sections IIVI and IVl to present the details of the calculation of both the virtual and real parts 
of the NLO QCD corrections to gg — > tih. Section E] also includes a discussion of the tree 
level (q, q)g —>■ tih + (q, q) processes. In Section IVII we explicitly show the factorization 
of the initial state infrared singularities into the gluon distribution functions, and finally 
summarize our results for the NLO inclusive total cross section for pp — > tih at the LHC in 
Eqs. f|KH|) and (|KK|) - (fH3|) . Finally, numerical results for the total cross section are presented 
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in Section IVlII We collect most of the technical details, including a list of box and pentagon 
integrals, in a series of Appendices. 



II. THE CALCULATION: GENERAL SETUP 

The inclusive total cross section for pp — > tih at 0(a^) can be written as: 
&nlo(pp -> tih) = 

J2yTJ~- I dxidx2 [^( x i^)Fj{ x 2,v)d-% LO (x 1 ,x 2 ,ii) + (1^2)] , (1) 

ij v 

where J-f are the NLO parton distribution functions (PDFs) for parton i in a proton, defined 
at a generic factorization scale fMf — fj,, and a l ^ LO is the 0(a^) parton-level total cross section 
for incoming partons i and j, made of the channels qq,gg — > tt/i and (q,q)g — > tih(q,q), 
and renormalized at an arbitrary scale /x r which we also take to be \i r = \i. Throughout 
this paper we will always assume the factorization and renormalization scales to be equal, 
\i r = fij = [A, unless differently specified. The partonic center-of-mass energy squared, s, is 
given in terms of the hadronic center-of-mass energy squared, s H , by s = x\x 2 s H . At the 
LHC center-of-mass energy the cross section is dominated by the gg initial state, although 
the other contributions cannot be neglected and are included in this calculation. 
We write the NLO parton-level total cross section a^ LO (xi, x 2 , /i) as: 

ajjj) 
An 

a^ (x h x 2 ,fi) + 8a% LO (x 1 ,x 2 ,n) , (2) 



^Loi X ^ X 2^) = « S 2 W /lo^i. ^1 + ^/; 3 i0 (ii,i 2 ,/i) 



where a s (fi) is the strong coupling constant renormalized at the arbitrary scale /i r = fi, 
o % l {x\,x 2 ,ii) is the 0(a 2 ) Born cross section, and 5a% LO {x\, x 2 , fi) consists of the 0(a s ) 
corrections to the Born cross sections for gg, qq — > tih and of the tree level (q, q)g — > tih(q, q) 
processes, including the effects of mass factorization (see Section |VlJ) . 5a^ LO {xi, x 2 , /i) can 
be written as the sum of two terms: 

5a i J LO (x 1 ,x 2 ,fx) = J d{PS 3 )Y^\A,irt{ij -> tih)\ 2 + J d{PS 4 )J2\A-eai{ij - tih + l)\ 2 

= K J irt( X l, X 2i + Klal( X li X 2, /i) , (3) 

where \ A V i r t(ij —> tih)\ 2 and \ A rea i(ij — > tih + l)\ 2 (for ij = qq, gg and l = g, or ij = qg, qg and 
l = q,q) are respectively the 0(a^,) terms of the squared matrix elements for the ij — > tih and 



ij — > tih + I processes, and ^2 indicates that they have been averaged over the initial state 
degrees of freedom and summed over the final state ones. Moreover, d(PS3) and d^PS^) in 
Eq. © denote the integration over the corresponding three and four-particle phase spaces 
respectively. The first term in Eq. (JHJ) represents the contribution of the virtual one gluon 
corrections to qq — > tih and gg — > tih, while the second one is due to the real one gluon and 
real one quark/antiquark emission, i.e. qq, gg — > tih + g and qg(qg) — > tih + q(q). 

The 0(a s ) virtual and real corrections to qq — > tih have been discussed in detail in 



Ref. 



21 1, and will not be repeated here. In the following sections we present the general 



structure of the 0(a s ) virtual and real corrections to gg — > tih. The contribution of the 
(q, q)g initiated process will be considered in Section^ when dealing with the real part of the 
0(a1) cross section. The results presented in the following sections have been obtained by 
two completely independent calculations, based on a combination of FORM and Maple 
codes in one case, and on the Mathematica based code Tracer [3^ in the other. The matrix 
elements squared for the tree level processes gg — > tih, gg — > tih+g, and (q, q)g — ► tih+(q, q) 
have been checked with Madgraph [31]. The numerical results presented in Section IV III have 
been obtained with two independent Fortran codes. 

Finally, we observe that the scale dependence of the total cross section at NLO is dictated 
by renormalization group arguments, and fN LO (xi, x 2, A 4 ) m Eq. © must be of the form: 

s 

with f[ 3 (xi,X2) given by: 

fi(x 1 ,x 2 ) = 2 < 4irb fl J o (x 1 ,x 2 ) - / dz 1 P lk (z 1 )f^ (x 1 z 1 ,x 2 ) 
{ k \- J p 

J dz 2 P jk (z 2 )f i I k (x 1 ,X2Z2) I , (5) 



+ 



where p=(2m t + Mh) 2 /s, P%j{z) denotes the lowest-order regulated Altarelli-Parisi splitting 



function 



32| of parton % into parton j, when j carries a fraction z of the momentum of 



parton i, (see e.g. Section |VJ), and 60 is determined by the one-loop renormalization group 
evolution of the strong coupling constant a s : 

(hAt ' ] -b a* + O(a*) , b = ±( 1 4N- 2 -n lf ) , (6) 



d\n(jj 2 ) u 3 v 81 ' u 4vr V 3 3 
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with iV = 3, the number of colors, and nj/ = 5, the number of light flavors. The origin of the 
terms in Eq. (JSJ will become manifest in Sections IIVI IVl and IVII when we describe in detail 
the calculation of both virtual and real 0(a s ) corrections. 

III. THE TREE LEVEL CROSS SECTION FOR gg -»■ tih 

The tree level amplitude for the process 

g A {qi) + g B {q2) -> t(p t ) + i(p' t ) + KPh) , 

where q\ + q 2 = Pt + v't + Ph and A, B denote the color of the incoming gluons, is obtained 
from the three classes of Feynman diagrams represented in Fig. Q identified as s— channel, 
t— channel, and u— channel diagrams respectively. We find it convenient to organize the 
color structure of both the tree level amplitude and the one-loop virtual amplitude in terms 
of only two color factors, one symmetric and one antisymmetric in the color indices of the 
initial gluons. Following this prescription, the tree level amplitude for gg — > tih can be 
written as: 

A = A% ab [T A ,T B }+A a b {T A ,T B } , (7) 

where T A > B = X A ' B /2 in terms of the Gell-Mann matrices A A,B Q- A$ and A™* correspond 
to the terms in the amplitude that are proportional respectively to the abelian (or symmetric) 
and non-abelian (or antisymmetric) color factors and are explicitly given by: 

K = \(Ao,t + Ao, u ) , A n ab = Ao >s + \(A ,t - A , u ) , (8) 

where Ao, s , Ao,t, and Ao, u are the amplitudes corresponding to the sum of the s— channel, 
t— channel, and u— channel tree level diagrams in Fig. [TJ ^4o,s, Ao t t, and Aq, u are given 
explicitly in Appendix 1X1 

Due to the orthogonality between symmetric and antisymmetric color factors, the tree 
level amplitude squared takes the very simple form: 

from which we can derive the LO partonic cross section, upon integration over the final state 
phase space: 

a % i {x u x 2 ,n) = J d(PS 3 )^2\A \ 2 (xi,x 2 ,ti) , (10) 

6 
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(N 2 - 1) (\A™ b \ 2 + \A* b \ 2 ) - — (N 2 - 1)\A\ 



(9) 
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FIG. 1: Feynman diagrams contributing to the tree level process gg — ► ti/i. The circled crosses 
indicate all possible insertions of the final state Higgs boson leg, each insertion corresponding to a 
different diagram. 



where the dependence of |^4o| 2 on X\ and x 2 (through s = x\x 2 sh) and on the renormalization 
scale \x (through a s (p))) has been made explicit. 

When averaging over the polarization states of the initial gluons, the polarization sum of 
the gluon polarization vectors, e^qi, Ai) and e„(q 2 , A2), has to be performed in such a way 
that only the physical (transverse) polarization states of the gluons contribute to the matrix 
element squared. We adopt the general prescription: 

where i = l,2 and the arbitrary vectors n« have to satisfy the relations: 

n i* E A i) e ^(^' A ») = > n i E A i)<(<?i, \) = ' (I 2 ) 

Ai=l,2 Ai=l,2 

together with nf 7^0 and n 1 ^n 2 We choose n\ = q 2 and n 2 = qi, such that: 

E A ») = + 2 M ^ 2A • (13) 

A i= l,2 

Finally, the entire calculation is performed using Feynman gauge for both internal and 
external gluons. 



IV. NLO VIRTUAL QCD CORRECTIONS TO gg -> tth: THE a 9 v ? rt CROSS SEC- 
TION. 

The 0(a s ) virtual corrections to the gg — > tth tree level process consist of the self-energy, 
vertex, box, and pentagon diagrams illustrated in Figs. |2E The 0{pP s ) contribution to the 
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FIG. 2: C(a s ) virtual corrections to 55 — ► it/i: self-energy diagrams. The shaded blobs denote 
standard one-loop QCD corrections to the gluon and top quark propagators respectively. The 
circled crosses denote all possible insertions of the final state Higgs boson leg, each insertion 
corresponding to a different diagram. All i-channel diagrams (labeled as s!f)) have corresponding 



u-channel diagrams. 



virtual amplitude squared of Eq. J3J) can then be written as: 

^2\A virt (gg - tih)\ 2 = J2Y1 { AoA ki + AoM,^ = EE 27 ^ { A ° A k,) > ( 14 ) 

where Aq is the tree level amplitude given in Eq. (|7J). while Ad^. denotes the amplitude for 
a class of virtual diagrams that only differ by the insertion of the final state Higgs boson 
leg, i.e. D^j = Ylk^id W1 th Di = Si,Vi, Bi, Pi, j = s,t,u, and k running over all possible 
Higgs boson insertions, as illustrated in Figs. EK3 

The amplitude of each virtual diagram {AD i3 ) is calculated as a linear combination of 
fundamental Dirac structures with coefficients that depend on both tensor and scalar one- 
loop Feynman integrals with up to five denominators. The tensor i nteg rals are further 
reduced in terms of scalar one-loop integrals using standard techniques 0,Q. The 0(a s ) 
virtual corrections to gg — > tth involve pentagon tensor integrals of rank higher than one, i.e. 
Feynman integrals with five denominators and more than one Lorentz tensor index. These 
pentagon tensor integrals are not present in the corresponding corrections for qq — > tih. 
This introduces a new difficulty in the calculation, due to the numerical instabilities that 

8 



o> 



X 



y(l,2) 
V 1,8 





v d,2) 
V 3,s 



8 nnnnnnnr9nr(p<^ . 




v (l,2) 
V 4,t 




4f 



7(U) 

'6,s 



t 




g mnnRnRP — ® — t 



v (l,2,3) 
V 7,t 



innnnnnn 



innnn 




v (l,2,3) 
V 8,t 



g nnnnnnnr — ® — 



v (l,2,3) 
V 9,t 



8 mnnr^ 



g mnnnnnr 1 



v (l,2) 
V io,t 



TtHRRRRH 



7RRR 




7(3,4) 
' 10,1 



FIG. 3: 0(a s ) virtual corrections to gg — ► ti/i: vertex diagrams. The shaded blobs denote standard 
one-loop QCD corrections to the ggg, gti, or /lit vertices respectively. The circled crosses denote all 
possible insertions of the final Higgs boson leg, each insertion corresponding to a different diagram. 
Diagrams with a closed fermion loop have to be counted twice, once for each orientation of the loop 
fermion line. All t-channel diagrams (labeled as have corresponding u-channel diagrams. 
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FIG. 4: C(a s ) virtual corrections to 55 — ► tth: box diagrams. The circled crosses denote all 
possible insertions of the final Higgs boson leg, each insertion corresponding to a different diagram. 
Diagrams with a closed fermion loop have to be counted twice, once for each orientation of the loop 
fermion line. All t-channel diagrams (labeled as B^)) have corresponding u-channel diagrams. 



may arise as a consequence of the proportionality of the tensor integral coefficients to higher 
powers of the inverse Gram determinant (GD]_pf the full gg — > tth phase space. Indeed, 
the standard techniques introduced in Refs. |33|, |3^ allow us to rewrite a tensor integral 
as a linear combination of the linearly independent tensor structures that can be built, 
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FIG. 5: 0(a s ) virtual corrections to — > tt/i: pentagon diagrams. The circled crosses denote all 
possible insertions of the final Higgs boson leg, each insertion corresponding to a different diagram. 
All i-channel diagrams (labeled as Pa) have corresponding u-channel diagrams. 



for a given tensor rank, out of the independent external momenta and the metric tensor. 
The coefficients of the linearly independent tensor structures can be found by solving a 
system of linear equations, one for each independent tensor structure. As a result, they 
are proportional to inverse powers of the so called Gram determinant (GD), of the form 
GD = det(pi-pj) with pi and pj generic independent external momenta (for i, j — 1, . . . ,4, 
since only four out of the five external momenta are independent). The higher the rank of the 
original tensor integral, the higher the inverse power of GD that appears in the coefficients 
of its tensor decomposition. 

To briefly illustrate the problem, we parameterize the Gram determinant in terms of the 
tth phase space variables as 

GD = - [S - (2 ^ + Mh)2] [Ml + (s - su) 2 - 2M 2 h (s + su)} s ^-sin 2 ^-sin 2 tf sin 2 9 , 

(15) 

where s = X\X2S H is the partonic center-of-mass energy squared, and the tth phase space 
has been expressed in terms of a time-like invariant s t i= (pt + p't) 2 , polar angles (6, 6 t t) and 
azimuthal angles (0, <f) t t) in the center-of-mass frames of the incoming gluons and of the tt 
pair, respectively. As can be seen in Eq. (j!5j). the Gram determinant vanishes when two 
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momenta become degenerate, i.e. at the boundaries of phase space. Near the boundary of 
phase space it can become arbitrary small, giving rise to spurious divergences which cause 
serious numerical difficulties, since they appear in various parts of the calculation that are 
normally numerically, not analytically, combined. In the case of a 2 — ► 3 process, this 
problem arises for pentagon tensor integrals, when all the independent external momenta 
are involved, and it becomes more serious for higher rank tensor integrals. The probability 
that the Monte Carlo integration hits a point close to the boundary of phase space is not 
negligible and these points cannot just be discarded. 

We use two methods to overcome this problem and find agreement within the statistical 
uncertainty of the Monte Carlo phase space integration. In the first method, we impose 
kinematic cuts to avoid the phase space regions where the Gram determinant vanishes, and 
then extrapolate from the numerically safe to the numerically unsafe region using different 
algorithms. We have used extrapolations based on polynomial or trigonometric functions. 
We have also reproduced the analytic dependence of each pentagon diagram on the Gram 
determinant, tested it in the safe region of phase space, and used it to extrapolate to the 
unsafe region. A phase space point is kept only if the true and the extrapolated results come 
very close to each other, after repeated iterations. Each extrapolation has been repeated 
imposing cuts on different kinematic variables, until a stable answer, independent of the 
kinematic cuts, can be found. The details of the extrapolation procedure are very technical 
and we do not think they can be of interest to this discussion. In the second method, 
after having interfered the pentagon amplitudes with the Born matrix element, we eliminate 
all pentagon tensor integrals by simplifying scalar products of the loop momentum in the 
numerator against the propagators in the denominator wherever possible. The resulting 
expressions are very large, but numerically very stable, and we have used them to confirm 
the results obtained using the extrapolation methods explained above. 

After the tensor integral reduction is performed, the fundamental building blocks are 
one-loop scalar integrals with up to five denominators. They may be finite or contain both 
ultraviolet (UV) and infrared (IR) divergences. The finite scalar integrals are evaluated using 
the method described in Ref. J2J] and cross checked with the numerical package FF Q|. 
The singular scalar integrals are calculated analytically by using dimensional regularization 
in d = 4 — 2e dimensions. The most difficult integrals arise from IR divergent pentagon 
diagrams with several external and internal massive particles. We calculate them as linear 
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combination of box integrals using the method of Ref. [2^, |23j and of Ref. |24|. Details of 
the box and pentagon scalar integrals used in this calculation are given in Appendix [Bj All 
other scalar integrals, with two or three denominators, are commonly found in the literature. 

Inserting all diagram contributions into Eq. (|14|) . we obtain the complete 0(a^) contri- 
bution to the virtual amplitude squared, and integrating over the final state phase space we 
calculate <5ff rt in Eq. Q. The UV singularities of the virtual cross section are regularized 
in d = 4 — 2e uv dimensions and renormalized by introducing a suitable set of counterterms, 
while the residual renormalization scale dependence is checked from first principles using 
renormalization group arguments. The detailed renormalization procedure adopted in this 
calculation is explained in Section II V Al The IR singularities of the virtual cross section 
are extracted in d = 4 — 2e IR dimensions and are cancelled by analogous singularities in the 
0(al) real cross section. The structure of the IR singular part of the virtual cross section is 
presented in Section IIVB1 while the IR singularities of the real cross section are discussed 
in Section The explicit cancellation of IR singularities in the total inclusive NLO cross 
section for gg — > tih is outlined in Sections El and IVII 

Finally, we note that the tree level amplitude Aq in Eq. ()14)) has generically to be con- 
sidered as the (i-dimensional tree level amplitude. This matters when the Ad { i amplitudes 
in Eq. (jl4J) are UV or IR divergent. Actually, as will be shown in the following, both UV 
and IR divergences are always proportional to the tree level amplitude or parts of it and 
they can be formally cancelled without having to explicitly specify the dimensionality of the 
tree level amplitude(s). After UV and IR singularities have been cancelled, everything is 
calculated in d = 4 dimensions. 



A. Virtual corrections: UV singularities and counterterms 

Self-energy and vertex one loop corrections to the tree level gg — > tih process give rise 
to UV divergences. These singularities are cancelled by a set of counterterms fixed by 
well defined renormalization conditions. As required by renormalization group arguments, 
the renormalization of the fundamental propagators and interaction vertices of the theory 
reduces to introducing counterterms for the external field wave functions of top quarks and 
gluons (SZ^ , 5Z 3 ), for the top mass (Sm t ), and for the strong coupling constant (5Z as ). The 
counterterm for the top quark Yukawa coupling, gtth = m t/v, coincides with the counterterm 
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for the top mass, since the SM Higgs vacuum expectation value v is not renormalized at one 
loop in QCD. 

By carefully grouping subsets of self-energy and vertex diagrams, we can factor out the UV 
singularities of the O(a^) virtual amplitude and write them in terms of the tree level partial 
amplitudes Ao tS , Aoj, and Ao tU introduced in Eq. (JSJ) and defined in Appendix^ According 
to the notation introduced in Figs. |2E we denote by D it j (with D = S, V, 2 = 1,2,..., and 
j = s, t, u) a class of diagrams with a given self-energy or vertex correction insertion, summed 
over all possible insertions of the external Higgs field, one for each different diagram. We 
now define Auy(AD i:j ) to be the UV pole part of the corresponding amplitude. Using this 
notation, we find 
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^«(- 

47T I 2 



A uv (Av 3)a + Av 9it + Av % J 
&uv(A S2 , a + A Sa>t + A Sa>u + A Si , t + -4 54 



1 

2iV 
a 



4 



A 



o , 
1 

2iV 



x 



(Ao + \(Ao !t - A ,u)[T A , T B ] + i(A,t + A 0tU ){T A , T B }^j , 



(16) 



where nif — 5 corresponds to the number of light quark flavors, N = 3 is the number of colors, 
M s and Mt are defined as: 

'47T// 2 ^ e 



+ e 



(17) 



and we have already included in the top quark self-energy diagrams the top mass countert- 
erm. 

We notice that some of the UV divergent virtual corrections (V% !S , Vj^m), and Vs, (*,«)), as 
well as and 5Z 3 in Eqs. (fTHj) and (fTUjl below, have also IR singularities. In this section 



14 



we limit the discussion to the UV singularities only, while the IR structure of these terms 
will be considered in Section IIVBI To this purpose we have explicitly denoted by e uv the 
pole parameter. 

The corresponding counterterms are defined as follows. For the external fields, we fix the 
wave-function renormalization constants of the external top quark fields using the on-shell 
subtraction scheme: 



(6Z?) 



uv 



Air 



N 
~2~ 



1 

2iV 



1 



4 



(18) 



while we renormalize the wave-function of external gluons in the MS subtraction scheme: 



(SZ- 



3)UV 



Air 



(4tt)T(1 + e) 



-N n 



if 



1 



1 



In 



11)7 



(19) 



.3 3 

according to which we also need to consider the insertion of a finite self-energy correction 
on the external gluon legs. This amounts to an extra contribution 



^ = £<*)T<1 + <)(§*-§«,,) in g) 



(20) 



which is important in order to obtain the correct scale dependence of the NLO cross section. 

We define the subtraction condition for the top-quark mass mt in such a way that m t is 
the pole mass, in which case the top-mass counterterm is given by: 



m t 



at* 



AlT 



N 



1 



-M t-^ — + 4 



2N 



(21) 



This counterterm has to be used twice: to renormalize the top-quark mass, in all diagrams 
that contain a top quark self-energy insertion, and to renormalize the top quark Yukawa 
coupling. As previously noted, the expressions in Eq. (jTSjl already include the top-mass 
counterterm. 



Finally for the renormalization of a s we use the MS scheme, modified to decouple the 
top quark 
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The first light flavors are subtracted using the MS scheme, while the 



divergences associated with the top-quark loop are subtracted at zero momentum: 



5Z as 



Air 



(4tt)T(1 + e) 



n lf -—N 



1 2 



1 1 / ^ 

h In ' ' 



mi 



(22) 



,3 1 3 

such that, in this scheme, the renormalized strong coupling constant a s (fi) evolves with 
riif = 5 light flavors. 
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Using the results in Eqs. ((TB)) - (J221) it is easy to verify that the UV pole part of <3ff rt : 



{Vvirt)UV-pole 



SlTlt 



(6ZV) uv + (8Z 3 ) vv + 8 uv + —± + 6Z 



(23) 



is free of UV singularities and has a residual renormalization scale dependence of the form: 



99 



-n lf + —N ] In ' 



2 , v 3 3-v-v. ' ' (24) 

as expected by renormalization group arguments (see the first term of Eq. (J3J). We note that 
the presence of s in the argument of the logarithm of Eq. (j^lj) has no particular relevance. 
Choosing a different argument would amount to reabsorbing some /i-independent logarithms 
in fl j of Eq. (|1J). 



B. Virtual corrections: IR singularities 

The structure of the IR singularities originating from the 0(a s ) virtual corrections to the 
tree level amplitude for gg — > tth is more involved than for the UV singularities. However 
it simplifies considerably when given at the level of the amplitude squared, and this is what 
we present in this section. 

The IR divergent part of the 0{pP s ) virtual amplitude squared of Eq. (|T^|) can be written 
in the following compact form: 

YY? ne **MdJ) = ( C ^ M vl + C ^vl + C ^vl) > (25) 

where Af t is defined in Eq. (fTTj) and we denote by A IR (ADij) the IR pole part of the amplitude 
of a given Djj class of diagrams. The result is organized in terms of leading and sub-leading 
color factors: 

N 2 

C X = ^(iV 2 - 1) , 

C 3 =(l + ^) (iV 2 - 1) , (26) 
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and the corresponding matrix elements squared M^), My e , and My\ are given by: 



M 



M 



M 



(i) 

V.c 



(2) 
V,e 



(3) 
V.c 



4 2 
-^r + — (-2 + A.) 



-IR C IR 



(\A% ab \ 2 + \A% b \ 2 ) 



1 

e IR 



+ — [(A n + A r2 ) | Ao, s + A , t \ 2 + (Ar 3 + A r4 ) \A ,u - A , s \ 2 ] , 



+ 



e 7fl e lR 

2 s t t-2mf 
1 s t t - 2m| b 2 



;-2 + A T1 + A T2 + a T3 +a T4 ; 



a£>|2 



A^-fKT + Kf) 



(27) 



where the IR nature of the pole terms has been made explicit. A^ b and AQ ab are defined in 
Eq. ©, while Ao <s , Ao jt , and Ao tU are given explicitly in Appendix |XJ Moreover, we have 
defined: 

'i + /V 



l\2 



Pti 



4m t A 1 

, Atf = In 



V 1 - %, 

and we have introduced the notation: A CT = \n(a/mf) and A Ti = Infa/ mf) where 



(28) 



a 


= (qi + Q2) 2 


J 






T\ 


= ™t - (<?i - 




= 2qi 


■Pt 


r 2 


= ™ 2 t ~ (l2 ~ 


Pt) 2 


= 2q 2 


■Pt 


T 3 


= ™t ~ (?a - 


P t ) 2 


= 2q 2 


■Pt 


T 4 


= - (<?i - 


Pt? 


= 2qi 


Pt 



(29) 



When we add the IR singularities coming from the counterterms that we have introduced 
in Section HV Al we can write the complete pole part of the IR singular 0{ot 3 s ) virtual cross 
section as: 



(±.99 \ 

\ u virt)IR-pole 



d(ps 3 ) J2 2ne **MdJ) + ^io {^ z 2 ] ) IR + w 



d{PS,)^-M^ (CiM% + C 2 M {2 \ + C 3 M 



®s tr ( 2 



-N + 



1 

N ) r, 



-a 



>i>! 
LO 



(30) 



As will be demonstrated in Section El the IR singularities of a 9 v f rt are cancelled by the 
corresponding IR singularities of offaz- 
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FIG. 6: Examples of C(a s ) real corrections to gg — > tih (first two diagrams) and of the tree level 
(QiQ)9 ~^ tih(q,q) processes (third diagram). The circled crosses denote all possible insertions of 
an external Higgs boson leg, each insertion corresponding to a different diagram. 

V. NLO REAL QCD CORRECTIONS TO gg tih: THE &™ al AND &™ al CROSS 
SECTIONS 

The NLO real cross section a 99 al in Eq. (j3J) corresponds to the 0(a s ) corrections to 
gg — > tih due to the emission of a real gluon, i.e. to the process gg — > tih + g, examples 
of which are illustrated in Fig. |U1 It contains IR singularities which cancel the analogous 
singularities present in the 0(a s ) virtual corrections (see Section II V Bj) and in the NLO 
parton distribution functions. These singularities can be either soft, when the energy of 
the emitted gluon becomes very small, or collinear, when the final state gluon is emitted 
collinear to one of the initial gluons. There is no collinear radiation from the final t and i 
quarks because they are massive. At the same order in a s , the of eai cross section corresponds 
to the tree level processes (q,q)g — > tih + (q,q), an example of which is also illustrated in 
Fig. El This part of the NLO cross section develops IR singularities entirely due to the 
collinear emission of a final state quark or antiquark from one of the initial state massless 
partons. The IR singularities can be conveniently isolated by slicing the gg — > tih + g and 
(q, q)g — > tih+(q, q) phase spaces into different regions defined by suitable cutoffs, a method 
which goes under the general name of Phase Space Slicing (PSS). The dependence on the 
arbitrary cutoff(s) introduced in slicing the phase space of the final state particles is not 
physical, and cancels at the level of the total real hadronic cross section, i.e. in a rea i, as well 
as at the level of the real cross section for each separate channel, i.e. in c 9 r 9 eal , o q9 eal , and o~l q eal . 
This cancellation constitutes an important check of the calculation and will be discussed in 
detail in Section fVTl 



18 



We have calculated the cross section for the processes 



g{qi) + gfo) -> t{ Pt ) + i( P ' t ) + h{ Ph ) + g(k) 



and 



(q, q)(qi) + 9{q%) -> t(p t ) + t(p' t ) + h(p h ) + (q, q)(k) 



with qi + q 2 = Pt + p' t + Ph + k, using two different implementations of the PSS method 
which we call the two-cutoff and one-cutoff methods respectively, depending on the number 
of cutoffs introduced. The two-cutoff implementation of the PSS method was originally 
developed to study QCD corrections to dihadron production 3^ and has since then been 
applied to a variety of processes (for a review see, e.g. 2^|). The one-cutoff PSS method 
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27| and extended to the case of massive quarks 



was dev elop ed for massless quarks in Ref. 
in Ref. \2^. 

In the next two sections we discuss the application of the PSS method to our case, 
using the two-cutoff implementation in Section IV Al and the one-cutoff implementation in 
Sect ion IVBl The results for a rea \ obtained using PSS with one or two cutoffs agree within the 
statistical errors of the Monte Carlo integration. In spite of the fact that both methods are 
realizations of the general idea of phase space slicing, they have very different characteristics 
and finding agreement between the two represents an important check of our calculation. 



A. Phase Space Slicing method with two cutoffs 

The general implementation of the PSS method using two cutoffs proceeds in two steps. 
First, by introducing an arbitrary small soft cutoff S s , we separate the overall integration of 
the gg — > tih + g phase space into two regions according to whether the energy of the final 
state gluon (k° = E g ) is soft, i.e. E g < 5 s y/s/2, or hard, i.e. E g > 8 s ^/s/2. The partonic real 
cross section of Eq. (J3J) can then be written as: 

= + , (31) 

where &soft ^ s obtained by integrating over the soft region of the gluon phase space, and 
contains all the IR soft divergences of o~ 9 r 9 eal . To isolate the remaining collinear divergences 
from dfa rd , we further split the integration over the gluon phase space according to whether 
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the final state gluon is {o 9 h 9 ard/coll ) or is not (v 9 h 9 ard/non _ coll ) emitted within an angle 9 from 
the initial state gluons such that (1 — cos#) < 5 C , for an arbitrary small collinear cutoff S c : 

-99 _ -99 i fr99 (oo\ 

u hard ~ u hard/coll ' u hard/non-coll ' V 3 *) 

In the same way, we isolate the collinear divergences in the cross section for the (q, q)g 
initiated processes and write the corresponding cross section as: 

^real = ®coll + ^non-coil ■ (33) 

The hard non-collinear part of the real gg- initiated cross section, ^ffard/non-coiv anc ^ ^ ne non ~ 
collinear part of the (q, q) ^-initiated cross section, <T 9 f on _ coll , are finite and can be computed 
numerically. 

On the other hand, in the soft and collinear regions the integration over the phase space 
of the emitted gluon or quark can be performed analytically, thus allowing us to isolate and 
extract the IR divergences of o 99 eoX and o q9 eoX . More details on the calculation of o 99 Q ^ t and 
^ffard are gi ven i n Section IV A II and Section IV A 21 respectively. The calculation of o q9 eal is 
described in Section IV A 31 

1. Real gluon emission, gg — ► tth + g: soft region 

The soft region of the phase space for the gluon emission process 

g%i) + 9 B (q2) - t( Pt ) + t{p' t ) + h( Ph ) + g c (k) (34) 
is defined by demanding that the energy of the emitted gluon (k° = E g ) satisfies the condition 

E g < S s ^ (35) 

for an arbitrary small value of the soft cutoff 8 S . In the soft limit (E g — > 0), the amplitude 
for this process can be written as: 

A sof t{gg -> tih + g) = 



Pi- 


e* 


?r 


€* 


Pi 


■k 


Qi 


■k 


p't- 


e* 


q2- 


e* 


p'i 


■k 


12 


■k 


Qv 


e* 


<?2' 


e* 


Qi 


■k 


12 


■k 



- T T T [^--^t) (A* + A-) - T T T ( ^- - ^- ) (A , u - A , s ) 
+ T A T C T B ri±—- l±-\ { Ao, l A,j 4 T»T< T ' [ Ail .. 



•CrpBrpA | 


( Pf 


e* 


q2- 


€* 




VPt 


■k 


Q2 


■k 


iBrpArpC | 


(P'i' 


e* 


Qi- 


e* 




U 


■k 


Qi 


■k 


'BrpCrpA f 


( 92' 


■e* 


Qr 


•e* 




V?2 


■k 


Qi 


•A- 



(36) 
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where A, B, and C are the color indices of the external gluons, while e M (&;, A) (for A = l, 2) is 
the polarization vector of the emitted soft gluon. Moreover, in the soft region the gg — ► tih+g 
phase space factorizes as: 

d(PS 4 )(gg - tth + g) d(PS 3 )(gg - tth)d(PS g ) soft 



(37) 

where d^PS^) and d(PS 3 ) have been defined in Section [TT1 while d(PS g ) so f t denotes the 
integration over the phase space of the soft gluon. Since the contribution of the soft gluon 
is now completely factorized, we can perform the integration over d(PS g ) so ft analytically, 
using dimensional regularization in d = 4 — 2e to extract the soft poles that will have to cancel 
the corresponding singularities in Eqs. (pffijl and (J27|) . The integrals that we have used to 
perform the integration over the phase space of the soft gluon are collected in Appendix O 
After squaring the soft amplitude A so ft, summing over the polarization of the radiated 
soft gluon, and integrating over the soft gluon momentum, the pole part of the parton level 
soft cross section reads 



pole 



(vfoftUe = j d(PS 3 ) (J d{PS g ) soft Y)A soft {gg^tth + gf 

= jd(PS 3 )^^(cMZ + C 2 M^ + C 3 M^) , (38) 

where Ci, C2, and C3 are defined in Eq. (J2HJ), while Adg e , M-s t -> an< ^ represent the IR 
pole parts of the corresponding matrix elements squared, and can be written as: 

M% = -M { il- 2 -{l + A\n{5 s ))(\Ar h \ 2 + \A$\ 2 ) , 

.Mg = -^g-^in^i^f+^i^rr+Ki 2 ) , 

M% = -M%+\\At b \\ (39) 

Note that in this section we do not explicitly denote the IR poles as poles in e IR , since all 
singularities present in a 9 r 9 e '^ 9 are of IR origin. 

After adding the IR divergent part of the parton level virtual cross section of Eq. (jHOj) we 
obtain: 



±.gg — (±.99 \ 1 (±.99 \ _ s \r 

(T s+V = {<J soft )pole + {(T virt )lR-pole ~ ~^ JV t 
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-4N\n(5 s )-^(llN-2n lf ) 



(40) 



where we can see that the IR poles of the parton level virtual cross section are exactly 
canceled by the corresponding singularities in the parton level soft gluon emission cross 
section. The residual divergences will be canceled by the soft+virtual part of the PDF 
counterterm when convoluting with the gluon PDFs as will be demonstrated in Section IV11 
The finite contribution to the parton level soft cross section is finally given by 



(A- 99 \ 

\ u soft J f mite 



d(PS 3 ) I / d(PS g ) soft Y)Asoft{gg^tth + g)\- 



finite 



d(PS 3 )^-Mj2 {CiMf + C 2 Mf + C 3 M 



(41) 



where the finite parts of the Ai^g\ ■> anc ^ A4g matrix element squared are explicitly 
given by: 
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(42) 



We note that s t t, Ptt, and A# are defined in Eq. (J25j) . while the function F(pi,pf) can be 
found in Appendix O (Eq. (J£3J)). 
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2. Real gluon emission, gg — ► tth + (7; hard region 

The hard region of the final state gluon phase space is defined by requiring that the energy 
of the emitted gluon is above a given threshold. As we discussed earlier, this is expressed 
by the condition that 

E 9 > 5 s ^f , (43) 

for an arbitrary small soft cutoff S s , which automatically assures that d~f^ ard does not contain 
soft singularities. However, a hard gluon can still give rise to singularities when it is emitted 
at a small angle, i.e. collinear, to a massless incoming or outgoing parton. In order to 
isolate these divergences and compute them analytically, we further divide the hard region 
of the gg —>■ tth + g phase space into a hard/ collinear and a hard /non- collinear region, by 
introducing a second small collinear cutoff S c . The hard/non- collinear region is defined by 
the condition that both 

-f— - > S c and -f = > S c (44) 

are true. The contribution from the hard/non- collinear region, o~l ard ^ non _ coll , is finite and 
we compute it numerically by using standard Monte Carlo integration techniques. 

In the hard/ 'collinear region, one of the conditions in Eq. ()44j) is not satisfied and the hard 
gluon is emitted collinear to one of the incoming partons. In this region, the initial state 
parton i with momentum q$ is considered to split into a hard parton %' and a collinear gluon 
g, % — > i'g, with qy = zqi and k — (1 — z)q{. The matrix element squared for ij tth + g 
factorizes into the Born matrix element squared and the unregulated Altarelli-Parisi splitting 
function Pu> = Pf v + ePL for i — > i'g, i.e.: 



\Ar e ai{ij - tth + g)\ 2 co <^ (4vra s )^| Ao(i'j - tih)] 2 ^ 1 , (45) 



where Ff v and P[ v denote the coefficients of the 0(1) and 0(e) parts of Pa>, while s ig = 2qi-k. 
In the case of gg — > tth + g the initial partons are gluons and the unregulated splitting 
function in d dimensions is (Pg g = 0): 

P ii ,( z )=P gg ( z )=2N(-^ + — + z{l-z)) . (46) 
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Moreover, in the collinear limit, the ij — > tth + g phase space also factorizes as: 



d{PS A ){ij -> tth + g) 
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d{PS z )(i'j^tth) 



zd^k 



(27r)( d - 1 )2Eg 
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(47T) 
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(47) 

where di^PS^) and d(PS 3 ) have been defined in Section HTJ while the integration range for s ig 
in the collinear region is given in terms of the collinear cutoff, and we have defined s' —Iq^-qy 
The integral over the collinear gluon degrees of freedom can then be performed analytically, 
and this allows us to explicitly extract the collinear singularities of & 99 ard as: 
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(4£ 



where i, i', and j are all gluons. As usual, these initial state collinear divergences are absorbed 
into the gluon distribution functions as will be described in detail in Section IV11 



3. The tree level processes (q, q)g — > tth + (q, q) 



The extraction of the collinear singularities of o q9 eal is done in the same way as described in 
Section IV A 21 for the gg initial state. In the collinear region, cos 9 iq > 1 — S C , the initial state 
parton i with momentum q^ is considered to split into a hard parton il and a collinear quark q, 
i — ► i'q, with qy =zqi and k — (1— z)qi. The matrix element squared for ij —>■ tih+q factorizes 
into the unregulated Altarelli-Parisi splitting functions in d dimensions: Pa> = P^ + eP^, and 
the corresponding Born matrix elements squared. The ij — *> tth + q phase space factorizes 
into the i'j — ► tth phase space and the phase space of the collinear quark. As a result, after 
integrating over the phase space of the collinear quark, the collinear singularity of <J q9 eal can 
be extracted as: 
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(49) 
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The collinear radiation of an antiquark in qg — > tth + q is treated analogously. In the case 
of (q, q)g — > tth + (q, q) we have two possible splittings: (q, q) — > g(q, q) and g — > qq. The 
0(1) and 0(e) parts of the corresponding splitting functions are: 

p; q (z) = \ (z 2 + a - -) 2 ) , 

P'(z) = -z(l-z) , 



Pi(z) 



N 2 — 1 /l + (1 
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2N 



N 2 - 1 

= -~2ir z ■ ( 5 °) 

Again, these initial state collinear divergences are absorbed into the parton distribution 
functions as will be described in detail in Section IVII 



B. Phase Space Slicing method with one cutoff 

An alternative way of isolating both soft and collinear singularities is to divide the phase 
space for the radiated parton into only two regions, according to whether all partons can be 
resolved (the hard region) or not (the infrared or ir region). In the case of gg — > tih + g, the 
hard and ir regions are defined according to whether the final state gluon can be resolved. 
The emitted gluon is considered as not resolved, and therefore part of the ir cross section, 
when 

s ig = 2pi ■ k < s min (Pi = qi, q 2 , Pt, Pt) , (51) 

for an arbitrary small cutoff s m i n , with k the final state gluon momentum which becomes 
soft or collinear. In the case of (q, q)g — > tth + (q, q), the emitted (anti)quark is considered 
as not resolved, and therefore part of the ir cross section, when 

s ig = 2pi ■ k < s min (Pi = qi,q 2 ,PuP' t ) , (52) 

with k the final state (anti) quark momentum which becomes collinear. The partonic real 
cross sections can then be written as the sum of two terms (ij = gg, qg): 

KLl = °l + Ohard > ( 53 ) 

where a} J r includes the IR singularities, both soft and collinear, while o % h d ard is finite. Following 
the general idea of PSS, we calculate a % ? r analytically, while we evaluate <J % l ard numerically, 

25 



using standard Monte Carlo integration techniques. Both af r and o^ ard depend on the cutoff 
s min , but the hadronic real cross section, a rea i, is cutoff independent after mass factorization, 
as will be shown in Section fVTl 

In order to calculate <5f/ we apply the formalism developed in Refs. 21, 
follows. 



as 



(a) We consider the crossed process 

h( Ph ) - t(pt) + t{p' t ) + g A {q x ) + g B (q 2 ) + g C (k) , (54) 

obtained from gg — > tth + g by crossing all the initial state colored partons to the final 
state, while crossing the Higgs boson to the initial state. All colored partons are hence 
considered as final state partons. For a systematic extraction of the IR singularities 
within the one-cutoff method, we organize the amplitude for h —>■ ggti+g, A r ^ 9tt+g , in 
terms of six colored ordered amplitudes, Aijk, which are the coefficients of all possible 
permutations of the color matrices T A ,T B ,T C , i.e.: 

A h r Zr +9 = E Ank^T*!*. (55) 

i,j,k = A,B,C 

The explicit color ordered amplitudes have very lengthy expressions and we do not give 
them in this paper. Since they are tree level amplitudes, they can be easily obtained. 
In the following we will however discuss in detail their properties in both the soft and 
collinear regions of the phase space of the extra emitted gluon. In this respect, we 
note that decomposing A^J 9tt+g in terms of color ordered amplitudes Aijk allows us 
to write the partonic real cross section as: 

h^ggtt+g _ I j/jdq \V^I A h ~^9gtt+g\2 



a real 



J d(PS 5 )Y,\A h r 2r tI+9 \ 2 , (56) 
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where 



5> 



h^ggtt+g 1 2 
real 



i,j,k=A,B,C 



\ Z + \A cba + Abac 4" *<4 S ca| + I^cas + -4cba + -4bca| 4" 



1-4-ABC 4" *-4 B , 4C + A ,\c 



\Acab 4" Acba 4" -4.^03 

1+1-4 ABC CA | 



+ C 3 



i,j,fc=A,B,C 



(57) 



is the full real amplitude squared, including both leading and subleading color factors 
(see Eq. (J2HJ) for a definition of Ci, C2, and C3). Each sub-amplitude squared in 
Eq. (JSTj) has very definite factorization properties in the soft or collinear regions of the 
phase space of the extra emitted gluon. 

(b) Using the one-cutoff PSS method and the factorization properties of soft and collinear 
divergences of the various amplitudes squared in Eq. (JSTj) . we extract the IR singu- 
larities from & r ^i 9tt+9 in d = 4 — 2e dimensions. In the soft and collinear limits we 



obtain: 



h^ggtt+g soft 



\h^ggtt+g\2 
oft I 



"h^ggtt+g collinear * h- 
a real > ° coll 



a h-fggtt+g 1 2 
^coll 



(58) 
(59) 



T l+9 = j d(PSMPSg)softJ2\^o 
= l d(PSMPSg)col^2\A, 

where we denote by d(PS g ) so f t (d(PS g ) co u) the phase space of the gluon g c in the 
soft (collinear) limit. In both the soft and the collinear limits, the cross section for 
h — > ggti + g integrated over the phase space of the IR singular gluon has the form: 

^fSoti 9 = J d(PS,)^Yl 

{d [K s , c (t; 1, 2; t) [AH + A$\ 2 + K SjC (t; 2, 1; t) \J® - A^ 
+C 2 [2K s ,c(t;i) (\A a b > ic) \ 2 + \A n ab ' M '^ 
4 (K SjC (t;l;t) + K s>c (t;2;t))\A: 
+ C 3 K s , c (t;i)\Al b > {c) \ 2 } , (60) 



l( c )|2 



+ 

ab,(c) 1 2 
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where the tree level amplitudes for the crossed process h — > ggtt, denoted by Aq s , A Qt , 
and Aq^ u , as well as their linear combinations A^'^ and A™ b '^ c \ can be obtained from 
the corresponding amplitudes for gg — > tth given explicitly in Appendix El by flipping 
momenta and helicities of the crossed particles. The functions K are either evaluated 
in the soft (Ks) or in the collinear limit (Kc), and will be explicitly given in Eqs. (jBHJ) 
and (|7U|) . Moreover, we notice that the arguments of the i^s.c functions are indices 
i— l,2,t,i denoting the external partons g A (q±), g B (q2), t(p t ), and i(p' t ) respectively. 
The explicit forms of both the pole and finite parts of &^ 99tt+g and d'^ l " tt+9 are 
given in Sections IV B II and IV B 21 respectively 

(c) Finally, the IR singular contribution of/ of Eq. (|53|) consists of two terms: 



a 



99_±99 + & 99 

ir 1 crossing \" / 



ir it ' crossing 



As described in Section IVB3I , a 99 is obtained by crossing g A , g B to the initial state 
and h to the final state in the sum of &^ 99tt+g and d-^," tt+9 , while d' 99 ossing corrects 
for the difference between the collinear gluon radiation from initial and final state 



partons 



271 ] . The IR singularities of a 9 f rt of Section IIVBI are exactly canceled 



by the corresponding singularities in a 99 . On the other hand, d' 99 ossing still contains 
collinear divergences that will be canceled by the PDF counterterms when the parton 
cross section is convoluted with the gluon PDFs, as we will show in Section fVTl 

When calculating the cross section for qg — ► tih + q in the collinear limit using the procedure 
described above, the resulting IR singular cross section a^ 9 is simply given by the initial state 



qg splitting 
e.g. Ref. 



unctions of Eq. (joTIj) convoluted with the corresponding Born cross sections (see, 

zi) 

a s 1 / 47r/z 2 



a q9 

ir 



' in in 



dz(l — z) e x 







27iT(l-e) V« 

[PM a 9 L 9 (g(q v )g(q 2 ) -> tth) + P gq (z) afMQi)^) - tth)] , (62) 

where the prime identifies the parton that originates from the splitting of a similar or different 
parent parton. The cross section for qg — > tih+q in the collinear limit is obtained in complete 
analogy with Eq. (pp2l . 

Finally, the hard part of the parton level cross section, & t ^ ard (ij =gg, qg, qg), is finite and 
can be calculated numerically. In this respect we note that, in the one cutoff method, the 
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soft and collinear limits of the real cross section, and consequently d'^ ard , are more sensitive 
to the smallness of the cutoff. A cut on the full invariant masses Si g is more drastic than 
two separate cuts on either the energy or the angle of emission of the extra gluon (q or q), 
and can be felt even by terms in the amplitude squared that do not contain singularities. 
These effects are very small, but large enough to affect the results at the level of precision 
of our calculation. It is therefore crucial, in particular for &l 9 ard , to model the Monte Carlo 
integration for each term in Eqs. (j56j) - (JoT)) separately, and to enforce term by term only the 
cuts on the Si g invariants that are actually present in each term. 



1. Real gluon emission h — > ggtt + g: soft region 

We first consider the case of soft singularities, when, in the limit of E g — >0 (soft limit), one 
or more Si g < s m i n . In the soft limit, the h — > ggtt + g phase space, as well as the full parton 
level real amplitude squared, factorize the dependence on the degrees of freedom of the soft 
emitted gluon, as illustrated in Eq. (p)%|) . The soft part of the parton level cross section can 
be calculated analytically according to Eq. (jHOJ). The soft limit of the K functions, K s , is 
explicitly given by: 



7 jt > 



K s (t; t) = *2l J d(PS g ) soft [f ti (g) + h 3 {g) + f ft {g)\ = S ti + S l3 + S j 

Ks(t;t) = ^ J d(PS g ) soft f a (g) = S a , (63) 

where — 1,2 denote the two external hard gluons with momenta q\ and q%. For any pair 
of partons (a, b) excluding the soft gluon, the soft f a b(g) functions introduced in Eq. (jHSJ) 
are given by: 

s agSbg Sag S bg 

with Sij = 2pi ■ pj both for massless and massive partons, and the corresponding integrated 
soft functions S a b are consequently defined to be: 

g 2 N f 

Sab = -^~ I d{PS g ) soft (a,b,g)f ab {g) . (65) 



/-(9)S J^._^_!^, ( 64 ) 
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'n the one-cutoff PSS method, the explicit form of the soft gluon phase space is given by 
23: 



d(PSg) aoft (a,b,- 



(Any AM) 



[SagSbgSab ~ ^l^ag ~ ^a S fe S ] * ds ag ds 



16tt 2 T(l - e 

^mtn Sag~)@(Smin ^bg) 



bg 



X 



(66) 



with A = s 2 ab — 4m 2 m 2 , while the integration boundaries for s ag and Sb g vary according to 
whether a and b are massive or massless quarks. 

The explicit form of the integrated soft functions S a b is obtained by carrying out the 
integration in Eq. (|65|) and analytic expressions for the S a b are collected in Appendix [DJ 
Finally, using Eq. (JBTTj) . Eqs. ([ ^ -(JBlj )) . and the results in Appendix |Dl the pole part of the 
parton level soft cross section can be written as: 



\ u soft J pole 



'soft 

E [ci 

+2C 2 



a. 



d(PS 4 )-^Af t x 

Z7T 



4 2 2 8 
-= + - + -A a - -A 

rnrn 

e z e e e 



|,C ( T + \A r 



- ( i 



: I ( - + --^A mi „) \X 



ab,(c) 1 2 



+ c 3 



1 - 



gtf - 2m\ 

Sttfitt 



ofe,(c)|2 




(67) 



while the corresponding finite contribution is: 



(J 



h^ggtt+g\ 



'soft ) finite = I d{PS±)^Nt X 

E {cJ(A^ - 4A ff A min - 4A min + 8A^ 



7T 



)(K bi(c) r+K 



nab,(c) 1 2 



V 2 ^ S *l S t2 / 



2 2 

mi m 



■ ( - \A\ + A T3 + A T4 + ^ + ^ ) |< - 4 C) 12 



+ C 2 



-2A r 



s t t - 2m\ 

Sttfitt 



At 



St2 S t i 

2m\ 



%t 



(Ja+Jb) \^ C> \ 2 + \< 



+4 ( -4A mm + AAl m - \^ - \Al - \A\ - \a\ - \A\ 



+A T1 + A T2 + A T3 + A T4 + 



m 



2 2 2 \ 

mi mf mi \ 
+—+—+— 

S*l -5*2 -5*1 S t 2 / 
,2 



+ c 3 



Sttfttt 



+ — K 



o6,(c) 1 2 



, s 4 f — 2m 2 \ m£ 
^A min ( 1 A t t J H — j= \ J a + ^> 



a6,(c)|2 



(68) 
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where a, Tj, A CT , and A Ti are defined in Eq. (j2Uj) an d right before it, s t t, /%, and A# are 
defined in Eq. (|2*8Jl , A t j is given in Eq. (jD5|) , A min is: 

A min = In {^^j , (69) 

and the functions J a and Jf, are given in Eq. (|D4J) . 



2. iZeaZ gluon emission h — > ggtt + collinear region 



We now turn to the case of collinear singularities, which arise when one of the two final 
state gluons i (i — g A , g B ) and the hard extra gluon g (g — g) become collinear and cluster 
to form a new parton i' (also a gluon), i + g — > i' , with the collinear kinematics: = zqy 
and k = {1 — z)qv . In the collinear limit, the h — > ggtt + g phase space as well as the full 
parton level real amplitude squared factorize the dependence on the degrees of freedom of 
the collinear emitted gluon, as illustrated in Eq. ([59)1 . The collinear part of the parton level 
cross section can be calculated analytically according to Eq. (|5U|) . The collinear limit of the 
K functions, K c , is explicitly given by: 



Kc(t;i,j;t) 



d(PS g ) 

coll 



^9 S jygg-n 



2tt r(i 



2 lJtj 
I /4vW 1 
e 



- Jf 1 ■ 2n lf f« 



1 99-^9 



If -J 



+ 1 



99-*g 



) 



K c {t;i;t) 



+2n lf I q ^ g (0,0) 
Ng; 



d(PS, 



g I coll " 



'JT ' + "iff 



qq^9] 



2vr r i-e 



e 



'mm "mw- 



L gg-*g\ 



+ ni f lgg_>g(0,0) 



K c (t;t) = 



(70) 



where — 1,2 denote the two external hard gluons with momenta gi and t^- 
In the one-cutoff PSS method, the collinear gluon phase space is 



d(PS g ) coU (i,j, z) 



(47T) 



i6tt 2 T{T^7) s ^ dSi9 ^ 1 ~ ^ tdze ( Sm ™ - ■ ( 71 ) 



The collinear functions f l Jf l are proportional to the Altarelli-Parisi splitting function for 
ig — > i' and explicitly factorize the corresponding collinear pole, i.e. [26. 27]]: 



rig^i 
J a 



' ab 



2 



(72) 



31 



where both i and i' are gluons and therefore Pi 9 -+i> corresponds to the P gg (z) splitting 
function given in Eq. (|46|). The lower indices a and b have been used to specify the integration 
boundaries on z. In order to avoid double counting between soft and collinear regions of 
the phase space of g , it is crucial to impose that only one Si g at a time becomes singular, 
i.e. satisfies the condition Si g <s m i n . The advantage of having reorganized the amplitude in 
terms of color ordered amplitudes, as in Eq. (|55|h is that the fH~^ 1 collinear functions have 
a very definite structure: they are all proportional to (saiSigSgb)" 1 , for a, b= g A ,g B ,t,i, and 
the integration boundaries are then found by imposing that: 

Sai ZS a i> Smin * Z 7> Z\ 



s gb = (1 - z)si> b > s min — ► z<l-z 2 = l- • (73) 

Si' & 

The terms proportional to f q ^ 9 come from the fact that a pair of collinear final state 
massless quarks (riif = 5 is the number of massless flavors) can also mimic a gluon. The 
corresponding collinear function is: 

fqq^g _ 2 Pgg-g{ z ) (7A \ 
1 N s qq ' [<V 

where both the 0(1) and (9(e) parts of the splitting function P qq -^ g are defined in Eq. (|5U|). 
Note that we do not attach any lower index to f q ^ 9 because the integration over z has no 
singularities and can be performed over the entire range from z = to z — 1. 



The analytic form of the integrated collinear functions I gg ^ g (zi, z 2 ) and I aa 



(0,0) is 



obtained by carrying out the integration in Eq. (|7U|) . and is explicitly given by 27]: 

W*i, *>.) = 7 + ^ - 2) - y + (y " ^) 6 + 0(e 2 ) , 
W0,0) = iQ + ^+OCe 2 ) . (75) 
Finally, using Eq. (fHUJ) and Eqs. (J70*j) - (J7H|) . the pole part of the parton level collinear cross 
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section can be written as: 



cr 



h^ggtt+g\ 
coll ) 



pole = / d(PSt)-±M t - X 



2tt e 
-4A ff + 8A min + 2( 



Aab,(c),2 , I jna6,(c)|2 



(An + A T2 ) 14? + - (A.3 + A T4 ) \AZ - ^ 
-4 (A T1 + A T2 + A T3 + A T J + 16A min + 4 f — - - 



|( c )|2 



t Cc) 1 2 



+ c 2 



while the corresponding finite contribution is: 



h-*ggtt+g^ 



coll 



I finite 



o:. s 



d(PS 4 )-^M 



2tt 



e a 



-2A* - 12AL„ + 8A CT A 



2A r 



11 _ 2r^ 

3 3 iV 



O i' 2 ^ _ 67 , 10n '/ 1 \ (\ A ab,(c)\2 , I j nab. 



+ ( -1a^ - \a% + 2A mm (A n + A T2 ) ) + < 



(c)|2 



+c 2 



ab, (c) 1 2 




+ ( -2 A ? 3 - %K + 2A mm (A T3 + A T4 ) ) \A® - A^ 2 



8A mm (A T1 + A T2 + A T3 + A r4 ) - 24A^ m - 2 (A£ + A^ 2 + + A?J 



11 2n lf 



3 A" 



-4A mm ----^f -4 -7^-^ + -^f 



2 2 67 10 m/ 



9 9 N 



\A 



ab,(c) 1 2 




(76) 



(77) 



where a, r iy A CT , and A r . are defined in Eq. (j2^|) and right before it, while A min is defined in 
Eq. 



Ml 



3. IR Singular Gluon Emission: Complete Result for a\ 

Summing both soft and collinear contributions to the h — > ggti + g cross section of 
Sections EBUEEl and crossing the final state gluons g A , g B to the initial state and the 
Higgs boson to the final state (which flips both helicities and momenta of these particles), 
yields of/ of Eq. (JoTj) as 

%99 _ (~h^ggti+g , ~h^ggtt+g\ / 7 o\ 
a ir -\ a soft + a coll )crossed ■ l'°J 
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The IR pole part of of r is given by 



(n yy \ 
\ u ir Jpole 



/ / a h^ggtt+g\ LI"" 

\\P'aoft Jpole i \Pcoll ) pole) crossed 



i .h-*ggtt+g\ 



Z7T 



-~n lf + ~N - — ) ~b% , 
3 v 3 N ) e LO ' 



(79) 



where M. 



—A4y e (see Eq. (J27Jl ) and therefore (<rf/)poZe completely cancels the IR sin- 
gularities of the virtual cross section (^f rt )iR-poie m Eq- PHjl . The IR finite part of of/ is 
given by 



finite 



((*. 



h^ggtt+g^ 
soft I 



. / a h^ggtt+g 
finite T I0" co » 



) finite ) 



crossed 



(80) 



with (v s ^t,wi}~ 9 ) finite given in Eqs. (JSBI and (jZ7|). 

Finally, as described in Section IV Bl the partonic cross section for the IR singular real 
gluon radiation for the process — > tt/i + g using the one-cutoff PSS method, of/, is 
obtained from of/ by adding o 9 ^. osaing (see Eq. (jHU))- The cross section cr 9 J, ossing accounts for 
the difference between initial and final state collinear gluon radiation and contributes to the 
hadronic cross section as 

" x dz 



a 



mi 



0J J X2 * * 

: 

/4vr/i 2 ^ e 



(81) 



2tt \s 
11 ln (/ 



6 



3 N 

z 



r(i-e 

.2 

— e 



x 



1 

e 

7r" 67 5 n^/ 
3~ ~ 18 + 9^V 

(l-zf- e 



5(1 -z) 



m-^ ]+ + -^ +z{1 - z n\ ' (82) 

in terms of regularized plus functions (see Ref. [27] for the exact definition). As will be 
demonstrated in Section IVTl these remaining IR singularities will be absorbed into the gluon 
PDFs when including the effects of mass factorization. 



VI. THE TOTAL CROSS SECTION FOR pp -» tth AT NLO QCD 

The total inclusive hadronic cross section for pp — > tth is the sum of the contribution 
from the gg initial state, the qq initial state and the (q, q)g initial states 

<?nlo(pP -> tth) = a 99 LO {pp -> tth) + af LO {pp -> tth) + <J% 9 LO {pp -> tth) . (83) 
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As described in Section HH <J^ LO {pp — > tth) is obtained by convoluting the parton level NLO 
cross section a^ LO (pp — > tth) with the NLO PDFs J-f(x, f) (* — 9> <?)> thereby absorbing the 
remaining initial state singularities of <5cx^ iC , into the renormalized PDFs. In the following we 
demonstrate in detail how this cancellation works in the case of the gg and ( q, q) q initiated 



processes. The case of the qq initiated process is discussed in Section V of Ref . |2l| , where we 
presented in detail the contribution of the qq initial state to a NLO {pp — > tth). a q ^ LO {j>p — > tth) 
can be obtained from there with obvious modifications, and will not be repeated here. 

First the parton level cross section is convoluted with the bare parton distribution func- 
tions J-"f(x) and subsequently the T\{x) are replaced by the renormalized parton distribution 
functions, J-f(x,fif), defined in some subtraction scheme at a given factorization scale fif. 



Using the MS scheme, the scale-dependent NLO parton distribution functions are given 
in terms of the bare Ff(x) and the QCD NLO parton distribution function counterterms 
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271 ] as follows: 



(a) For the case where an initial state gluon, quark or antiquark {b = g,(q,q)) split 
respectively into a qq or (q, q)g pair {b' = (q, q),g): 



2tt 



ra- 



1 dz 

z 



(84) 



for both the one-cutoff and two-cutoff PSS methods, where is defined in Eq. (J3UJ). 
(b) For the case of g — > gg splitting: 



b.l) two-cutoff PSS method: 



&{x,p,) = TUx 




aiJ } r(i-c) 



l-Ss 



dz 

z 



where P gg is Altarelli-Parisi splitting function given in Eq. pBjl . 
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b.2) one-cutoff PSS method : 



+ 




2tt I (j? f J T(l -e) \ej \6 3 N 

1 dz 



d i (4) (f 



(86) 



where igg is the regulated Altarelli-Parisi splitting function given by: 



O^^^ + i^ + ;(!-,)) . (87) 

The 0(ck s ) terms in the previous equations are calculated from the 0(a s ) corrections to the 
b — > b'j splittings, in the PSS formalism. Moreover, note that in Eqs. (|84 |l -(|8f) )l we have 
carefully separated the dependence on the factorization (/ij) and renormalization scale (/x r ). 
It is understood that a s = a s (fi r ). The definition of the subtracted PDFs is indeed the only 
place where both scales play a role, and the only place where we have a dependence on [if. 
In the rest of this paper we have always set \x r = fif = // and we will also give the master 
formulas for the total NLO cross section, Eqs. (|8*%|) - (jn3j) . using /i r = /i/ = yU. We have checked 
that this simplifying assumption has a negligible effect on our results and we will comment 
more about this in Section fVIII 

In the two-cutoff PSS method, when convoluting the parton gg cross section with 
the renormalized gluon distribution function of Eq. (}85|) . the IR singular counterterm 
of Eq. exactly cancels the remaining IR poles of a 9 v f rt + of^ t and (^hard/coil- ^ ne 
one-cutoff PSS method, the IR singular counterterm of Eq. ((HE]) exactly cancels the IR 
poles of (^crossing- Finally, the complete 0(a^) inclusive total cross section for pp — > tth in 



the MS factorization scheme when only the gg initial state is included, i.e. <J% 9 LO (pp — > tth) 
of Eq. (jHB^) . can be written as follows: 
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1) two-cutoff PSS method 
1 



a 



99 



- J dx x dx 2 {j^(xi,/i)JJ(x2,Ai) [dj% } (x 1 ,x 2 ,n) + (a^ rt ) /inite (xi, x 2 , /i) 

+ (c r soft)finite(Xl } X 2 , fi) + a 9 s 9 +v+ct (x U X 2 , //) + (!<-> 2)] } 



+ - / dx\dx2 



l-<5« 



.t'l 



<iz 
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x CTf»(ari,x 2 ,/i)^ln f ^— — — - 



+ 



27T 



: 2 



P OT (z) + (1^2)} 



(7 



hard/non—coll 



(x 1 ,x 2 ,/i) + (l ^2)} 



(88) 



where of+^+ct is obtained from the sum of (<T^f rt )c/y-poie of Eq. (J2HJ) , of+t, of Eq. pUJ) . 
and the PDF counterterm in Eq. ()85|) as follows 



a 



99 

s+v-\-ct 



2^ 



11 2riif , „ ri / m? 

N ± + 4N\n(5 s ) In' * 
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(7 



99 



(89) 



2) one-cutoff PSS method 

= ^ / <ixi(ix 2 {j^(x 1 ,/j,)J^(x 2 ,fj l ) [af (x 1 ,x 2 ,iJ,) + {a^ rt ) f in u e {xi, x 2 , fJ.) 

/*) + (!<- 2)]} 
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1 dz 
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xa^ix^x^ii)— 2N\n 
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(1-z). 



z 

1 — z 
+ + z(l-z) 
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+ (1 <-> 2) 

+ U dXldX2 {L d i 

x^(x 1>a ; 2 ,/i)^2iV 



+ (1 <-> 2) 
1 
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1 - 2 
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+ ^(1 - z) ln(l - z) +z 



ln(l - z) 
1 - z 



+ 



cbiofe 2 {J^{x 1 ,fj l )J^(x 2 ,fJ>) af^ rd {x u x 2 , n) + (1^2)} 



(90) 



where a 99 +ir+ct is obtained from the sum of {a 9 v 9 irt )uv- P oie of Eq. ((23), \p" irt )m- po u 



<!9 



of 



Eq. (jSHJ), (o r f/) P oie of Eq. (fT9*|) . the part proportional to 5(1 — z) of a 
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9 J ossing ofEq. dHU), 



and the PDF counterterm in Eq. (|86|). and can be written as: 

_2 



a 



99 

v+ir+ct 



2^ 



HiV - 2n 

3 3 



+ 2iV 1 £ 

1 3 18 9N 



a 



99 
LO ■ 



(91) 



We note that aff LO is finite, since, after mass factorization, both soft and collinear singu- 
larities have been canceled between a 9 f rt + o 99 Q ^ t and d' 9 , 9 ard / cM in the two-cutoff PSS method, 
and between a 9 ? rt and a 99 in the one-cutoff PSS method. The last terms respectively de- 
scribe the finite real gluon emission of Eq. ()32|) and (|53|). Note that when collecting all the 
terms in Eqs. (J8*8)l and (J9*U|) that are proportional to ln(/i 2 /s), one obtains exactly the last 
two terms in Eq. (j^J), as predicted by renormalization group arguments. 

For the (q, q)g initiated processes we find 



1) two-cutoff PSS method 

^2 f dxxdxA I —J : f(—,fM)J^(x 2 ,iJ,)x 
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NLO 



i=q,q 

a 9 L 9 o(xi,x 2 ,n) 



z 



J xi z z 



s (1 - z) 2 5 f 



+ (1 - 2) 



+ / dx l dx ? {F?( X l^)Fg( X 2^)°non-cou( X li X 2^) + (1 2)} 

i=q,q 



(92) 



2) one-cutoff PSS method 



NLO 



2 1 X! / ^1^2 < / —J^{—,fi)Jf{x 2 ,fj,) 
dz ^1 / *^ i 



^(z)ln 



. (1 ~ f) 



j xi z z 

af {x l ,x 2 ,n) 



P^(z)\n 



:i-z) 



+ (1 <-> 2) 



+ Yl j dx l dx 2 {F?( X l,V)Fg( X 2,^)vZrd( X ^ X 2^) + i 1 2 )} 



l=q,q 



(93) 
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where Pfj and are the 0(1) and (9(e) contributions to the splitting functions as 
given in Eq. (joHj) . The last terms respectively describe the finite gluon/quark emission 
of Eqs. (O and (JS3J). 

We would like to conclude this section by showing explicitly that the total NLO cross 
section, a NLO , does not depend on the arbitrary cutoff's introduced by the PSS method, i.e. 
on s m i n for the one-cutoff method and on 5 S and S c for the two-cutoff method. The cancel- 
lation of the PSS cutoff dependence is realized in a rea i by matching contributions that are 
calculated either analytically, in the IR- unsafe region below the cutoff(s), or numerically, 
in the IR-safe region above the cutoff (s). While the analytical calculation in the IR-unsafe 
region reproduces the form of the cross section in the soft or collinear limits and is therefore 
only accurate for small values of the cutoff (s), the numerical integration in the IR-safe region 
becomes unstable for very small values of the cutoff(s). Therefore, obtaining a convincing 
cutoff independence involves a delicate balance between the previous antagonistic require- 
ments and ultimately dictates the choice of neither too large nor too small values for the 
cutoff(s). The Monte Carlo phase space integration has been performed using the adap- 
tive multi-dimensional integration program VEGAS j3] as well as multichannel integration 
techniques HHQ- 

In Figs. [Hand El we consider the two-cutoff PSS method and study the independence of 
cr NLO (pp — > tth) on 5 S and 5 C separately, by varying only one of the two cutoffs while the 
other is kept fixed. In Fig. [?[ 5 S is varied between 10~ 5 and 10~ 3 with 5 C = 10~ 5 , while in 
Fig- El $c is varied between 10 -6 and 10~ 4 with 5 S = 10~ 4 . In both plots, we show in the upper 
window the overall cutoff dependence cancellation between cr 9 s 9 oft +(T g h 9 ard/coll and a 9 h g ard/non _ coll 
in <y 99 al - We do not include the corresponding contributions from the Born and the virtual 
cross sections since they are, of course, cutoff independent. Similar plots could be obtained 
for the other two sub-channels, qq and qg + qg. We illustrate the point using just the gg 
channel, since the qq channel has already been thoroughly studied in Ref. |2l|, while the 
cutoff dependence of the qg + qg channel is trivial. In the lower window of the same plots 
we complement this information by reproducing the full <J NLO , including all channels, on a 
larger scale that magnifies the details of the cutoff dependence cancellation. The statistical 
errors from the Monte Carlo phase space integration are also shown. Both Figs. U\ and El 
show a clear plateau over a wide range of 6 S and 5 C and the NLO cross section is proven to 
be cutoff independent. The results presented in Section IVIII have been obtained by using 
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FIG. 7: Dependence of a NLO (pp — ► ii/i) on the soft cutoff <5 S of the two-cutoff PSS method, at 
ys^" = 14 TeV, for M h = 120 GeV, n = m t + M h /2, and 5 C = 10~ 5 . The upper plot shows the 
cancellation of the ^-dependence between o 99 ^ + a hard/coii an< ^ a hard/non-coir ^ ne l° wer pl°t 
shows, on an enlarged scale, the dependence of the full a NLO = o 99 LO + a qq LO + a q9 LO on 5 S with the 
corresponding statistical errors. 



the two-cutoff PSS method with 5 S = 10~ A and 5 C = 10~ 5 . 

We now turn to the one-cutoff PSS method and, following the same criterion adopted 
for the case of the two-cutoff PSS method, we summarize in the upper window of Fig. El the 
behavior of the different cutoff dependent contributions to the real gg — > tih cross section, 
i.e. a 99 and cr 9 ^ ard1 as well as the resulting cutoff independence of cr 99 al . The lower window 
of Fig. El shows the full cr NLO , where all tih subprocesses are included, on an enlarged scale. 
The statistical deviations due to the Monte Carlo integration are also shown, and therefore 
the stability of the integration procedure can be appreciated. In Fig. El s m i n is varied over 
several orders of magnitude and the presence of a clear plateau over most of the s m i n range 
is evident. The results presented in Section IVlIl have been cross-checked using the one-cutoff 
PSS method with s m „ = 10 GeV 2 . 
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FIG. 8: Dependence of a NLO (pp — > tth) on the collinear cutoff 5 C of the two-cutoff PSS method, 
at v /s^ = 14 TeV, for M fe = 120 GeV, fj, = m t + M h /2, and 6 S = 10~ 4 . The upper plot shows the 



cancellation of the ^-dependence between a 



99 
soft 



+ a 



99 

hard/coll 



, and a 



99 

hard /non— coll ' 



The lower plot 



shows, on an enlarged scale, the dependence of the full <j nlo = a 
corresponding statistical errors. 



nlo + <T q NLO + vVlo on S c with the 



VII. NUMERICAL RESULTS 

In this section we summarize the most important numerical results for a NLO {pp — > till) 
and illustrate the impact of NLO QCD corrections on the tree level cross section. In partic- 
ular, we discuss the renormalization/factorization scale dependence of a NLO with respect to 
a LO , and illustrate the dependence of both LO and NLO cross sections on the Higgs boson 
mass. Results for a LO are obtained using the 1-loop evolution of a s (^i) and CTEQ5L parton 
distribution functions Q|, while results for a NLO are obtained using the 2- loop evolution 
of a s (fi) and CTEQ5M parton distribution functions, with a^ LO {Mz) = 0.118. According 
to the renormalization prescription adopted in this paper and explained in Section II V Al 
throughout our calculation we use for the input parameter m t the top quark pole mass. 
Results are presented for mt — 174 GeV, while the uncertainty introduced by varying mt 
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FIG. 9: Dependence of a NLO (pp — > tth) on the s m i n cutoff of the one-cutoff PSS method, at 
yj~s~^ = 14 TeV, for = 120 GeV, and fj, = mt + The upper plot shows the cancellation 

of the s m j n -dependence between af® and cr 9 ^ ard . The lower plot shows, on an enlarged scale, the 
dependence of the full a NLO = a g N g LO + a qq LO + <J qq LO on s m i n with the corresponding statistical 
errors. 



within its experimental uncertainty is discussed later in this section. We define the top 
quark Yukawa coupling to be gtth = m t/v where v— (GfVz)^ 1 ^ 2 = 246 GeV is the vacuum 
expectation value of the SM Higgs field, given in terms of the Fermi constant Gp. 

In Fig. [T01 we illustrate the dependence of both a NLO and a LO on the renormalization and 
factorization scales when the two scales are identical, i.e. when We have also 

studied the behavior of o NLO when the renormalization and factorization scales are varied 
independently and noticed no appreciable difference with respect to the case in which the 
two scales are identical. This justifies our decision to present results only for 
We also illustrate in Fig. the \x dependence of the NLO cross section for each parton 
level channel independently. We use Mh = 120 GeV for the purpose of these plots. As 
expected, Fig. El shows that the NLO cross section has a much weaker scale dependence 
and represents a much more stable theoretical prediction. In Fig.[T21 we plot u LO (pp —> tth) 
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FIG. 10: Dependence of cr LOiNLO (pp — > tth) on the renormalization/factorization scale /i, at */s 
14 TeV, for Af fc = 120 GeV. 



and cr NLO (pp — > tt/i) as functions of the Higgs boson mass, for ^^1^= 14 TeV and two values 
of the common renormalization/factorization scale, fi = m t + and fi = 2m t + M^. We 

consider 100 GeV < < 200 GeV since the production of a Higgs boson in association 
with a pair of top quarks at the LHC will play a crucial role only for relatively light Higgs 
bosons. The information gathered from this plot nicely complements what has already been 
shown in Fig. We summarize a sample of results from both Figs. El and El i n Table |H 
where we also provide the LO cross section, W LO , calculated using the 2- loop evolution of 
a s (n) and CTEQ5M PDFs. This can be useful to separately evaluate the impact of the 
full set of NLO QCD corrections as opposed to the subset of them that mainly correspond 
to the NLO running of a s (/i). The error we quote on our values is the statistical error of 
the numerical integration involved in evaluating the total cross section. We estimate the 
remaining theoretical uncertainty on the NLO results to be of the order of 15 — 20%. This 
is mainly due to: the left over /i-dependence (about 15%), the dependence on the PDFs 
(about 6%), and the error on m t (about 7%) which particularly plays a role in the top quark 
Yukawa coupling. 

It can also be useful to quote the impact of NLO corrections in terms of a so called 
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FIG. 11: Dependence of cr NLO (gg, qq, qg + qg — > tth) on the renormalization/factorization scale ji, 
at J^=IA TeV, for M h = l2Q GeV. 



.fT-factor, defined as the ratio between the NLO and LO cross sections: 

OlcAw -> tth) 

We can see in Fig. El that, for a SM Higgs boson of mass Mh — 120 GeV, the i^-factor for 
pp — > tth is larger than unity when /i>0.4/i for fi = m t + Therefore, over a broad 

range of the commonly used renormalization/factorization scales, NLO QCD corrections 
increase the LO cross section. Using the results of Table HI the K-factors for a sample of 
Higgs boson masses and renormalization/factorization scales can easily be calculated, both 
using a LO and a L0 . We notice, however, that the fT-factor defined in Eq. (jMjl is affected by 
a very strong scale dependence, the same as a LO , Therefore, when the .fT-factor is used to 
obtain a NLO from must be used in matching a LO and K corresponding to the same 

/x-scale. 

VIII. CONCLUSIONS 



The associated production of a Higgs boson with a pair of top quarks will play a very 
important role in the discovery of a low mass Higgs boson at the LHC with a center of 
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FIG. 12: cr NLO (pp — > ti/t) and a LO (pp — ► it/i) as functions of M^, at y / i^~ = 14 TeV, for // = 
m f + M^/2 and fi = 2m t + M h . 

mass energy of = 14 TeV. With the statistics expected at the LHC, pp — » tt/i, with 
ft — > bb,r + T~ ,W + W~ ,77 will also play a crucial role in determining the couplings of the 
discovered Higgs boson, and will give the only handle on a direct measurement of the top 
quark Yukawa coupling. 

In this paper the inclusive cross section for pp — > tih production has been calculated, in 
the Standard Model, including full NLO QCD corrections. The NLO cross section shows 
a drastically reduced renormalization and factorization scale dependence, and leads to in- 
creased confidence in predictions based on these results. The overall uncertainty on the 
theoretical prediction, including the errors coming from parton distribution functions and 
the top quark mass, is reduced to only 15-20%, as opposed to the 100-200% uncertainty of 
the LO cross section. Including NLO QCD corrections increases the LO cross section for a 
broad range of commonly used renormalization and factorization scales, and over the entire 
Higgs boson mass range considered in this paper. This is summarized by saying that the 
f^-factor for renormalization and factorization scales in the range m t <fi< 2m t + Mh and 
Higgs boson masses in the range 100 GeV< Mh <200 GeV is between 1.2 and 1.6. 

The calculation of the NLO QCD cross section for pp — > tih contains several technical 
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M h (GeV) 


/" 


<J LO (fb) 


° LO (fb) 


o-jvio (fb) 


120 


mt 

m t + M h /2 

2m t 
2m t + M h 


582.92 ± 0.06 
520.47 ± 0.06 
450.09 ± 0.05 
405.54 ± 0.04 


616.81 ± 0.07 
553.25 ± 0.06 
480.80 ± 0.05 
434.59 ± 0.05 


718.64 ± 3.71 
697.27 ± 3.20 
662.66 ± 2.77 
634.36 ± 2.34 


150 


mt 

m t + M h /2 

2m t 
2m t + M h 


316.27 ± 0.03 
275.44 ± 0.03 
243.47 ± 0.03 
214.43 ± 0.02 


336.41 ± 0.04 
294.35 ± 0.03 
261.03 ± 0.03 
230.60 ± 0.02 


380.95 ± 1.81 
367.38 ± 1.52 
352.71 ± 1.35 
334.48 ± 1.18 


180 


mt 

m t + M h /2 

2m t 
2m t + M h 


187.44 ± 0.02 
159.32 ± 0.02 
143.77 ± 0.02 
123.85 ± 0.01 


200.46 ± 0.02 
171.15 ± 0.02 
154.74 ± 0.02 
133.65 ± 0.02 


221.63 ± 1.01 
214.01 ± 0.85 
206.59 ± 0.77 
194.42 ± 0.70 



TABLE I: Values of both cr LO (pp — > tth), a LO (pp — > tt/i), and cr NLO (pp — > tth), at - v /s^"=14 TeV, 
for a sample of different values of Mh and of the renormalization/factorization scales /i = /x r = / uy. 

difficulties that have been thoroughly explained in this paper (see also Ref. |2ll|). The 
NLO virtual corrections involve pentagon diagrams and consequently require the evaluation 
of both scalar and tensor pentagon integrals with several external and internal massive 
particles. Detailed information about the method used as well as explicit results for all the IR 
singular integrals appearing in the calculation are presented in a series of Appendices. Tensor 
pentagon integrals are affected by numerical instabilities and we discuss in this paper how we 
have calculated them in a numerically stable form. The NLO real corrections are complicated 
by the presence of IR divergences. We have calculated them in two different variations of 
the phase space slicing method, involving one or two arbitrary cutoffs respectively. The 
correspondence between the two phase space slicing methods is made explicit, and the 
agreement between them constitutes a powerful check of the technicalities used in their 
implementations. The techniques developed in this paper and in Ref. can now be 
applied to similar higher order calculations, in particular to the case of the associated bbh 
production at both the Tevatron and the LHC. 



46 



Acknowledgments 



We thank U. Baur, Z. Bern, and F. 
We are grateful to the authors of Ref. 



'aige for valuable discussions and encouragement. 



171 ] for a detailed comparison of the results. The 



work of S.D. (C.J., L.H.O., L.R.) is supported in part by the U.S. Department of Energy un- 
der grant DE-AC02-76CH00016 (DE-FG02-97ER41022, DE-FG-02-91ER40685, DE-FG02- 
97ER41022). 



APPENDIX A: TREE LEVEL AMPLITUDE FOR gg -► tth 

The amplitudes Aq^ Si Ao t t, and Ao, u introduced in Section [T] can be written as: 

A),s = W 2 S 9ith^(qi) £1/(92) u t A^ s Vt , 
A,t = ig 2 s g t the^qi)e v (q2)u t A^ t vt , 

Ao, u = ig 2 s 9tth^Mi)^(l2)u t A^ u vt , (Al) 

where g t ih = m t /v is the top quark Yukawa coupling, with v = 246 GeV the SM Higgs 
boson vacuum expectation value, while Aq" s , Aqj, and Aq" u represent the total s— channel, 
t— channel, and u— channel amplitudes, corresponding to the diagrams in Fig. More 
explicitly: 

/!/•"-' _ y((V 1 /I (2)|» I/ 1 /I (3)^ / A Q\ 
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where 



with 



*(i),i*> _ 1 A+ fa + m t ™ 

S [{Pt +PhF ~ m t\ 
A2),»v 1 ~ A~ A + 

A+ A + ra f M ^2- A + m f v 

°' f ^' [(Pt+Pfc) 2 -^] 7 [(fc-p0 2 -™?] 7 

7 " [(ft - ?i) 2 - m?] [(g 2 - pQ 2 - m 2 ] 7 ' 

,(3),^ _ A- k\ + m t „ - A- A + m t 



2 



7 1(pt-?i) 2 -m 2 ] 7 [(p£+P/i) 
AT = <^(/i -^9i-9 2 ) , (A3) 



^ Q = (?i - <? 2 ) Q <T + (gi + 2g 2 ) V° - (2gi + g 2 ) V a , 
are the individual amplitudes for the s— channel, £— channel, and u— channel diagrams in 

Fig.m 



APPENDIX B: BOX AND PENTAGON INTEGRALS 

We label the various one-loop box and pentagon scalar and tensor integrals appearing in 
the calculation of the 0(a s ) virtual corrections to 

g{qi) + gfa) -> t{ Pt ) + t(p' t ) + h{ Ph ) 

according to the diagram where they are encountered. Moreover we denote by DO, Dl p , 
D2 pv ', and D3 flUf> the scalar and tensor box integrals with one, two, and three tensor indices, 
and by EO, El 11 , E2 pv , and E?> pup the analogous scalar and tensor pentagon integrals. With 
this convention D0„(&), for instance, is the scalar box integral appearing in box diagram 
D^ k j , as labeled in Fig. 0] The external momenta are labeled as shown above, where qi,q2 
are incoming and pt,p' t , Ph are outgoing momenta with q\ + q 2 = Pt+Pt+Ph- It is convenient 
to express our results in terms of the kinematic invariants of Eq. (I29j) and: 

ui = (Pt+Ph) 2 - m 2 t , 

^2 = (p't+Ph) 2 -m 2 t . (Bl) 
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These kinematic invariants do not form a linearly independent set, but are related by: 



r 3 = a — Ti — 0J2 and r 4 = a — r 2 — u>\ . (B2) 

We also make frequent use of the shorthand notation A a = \n(a/m 2 ) with a = a, Ti,Ui. 

In the following, we explicitly give only the box and pentagon integrals that contain 
IR divergences. The IR divergences are extracted using dimensional regularization with 
d = 4 — 2e. We only give results for integrals arising from the s— channel and t— channel 
diagrams. The integrals for the u— channel diagrams can be obtained from the integrals of 
the corresponding t— channel diagrams by exchanging q\ «-> q 2 , i.e. by exchanging T\ «-> r 3 
and t 2 T4. The IR finite scalar integrals are evaluated by implementing the method 



described in Ref. 



24j and are cross checked against the FF package 



[35]. 



1. Box integrals 

(k) 

The scalar and tensor box integrals arising in the computation of box diagram D\ ■ are 
of the following form: 

W? • (B3) 

where 

N^ie-ml) , N 2 = (k + Pl ) 2 -ml , 
N 3 = {k + vi +P2) 2 -ml , N 4 = {k + Pl + p 2 + p 3 ) 2 - m 2 3 , (B4) 

Pi) P2, P3, and P4 = — pi — p 2 — P3 are the external (incoming) momenta connected to the 
box topology, and mo, mi, m 2 , and m 3 are the masses of the propagators in the box loop. 
We write the tensor integrals as a linear combination of the linearly independent tensor 
structures built of the independent external momenta p 2 , and p% plus the metric tensor 
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g pu . Our notation for the box tensor integrals is as follows: 

DV = D^pI + D< 2) p£ + D^pZ , 

D2T = D®gT + D^plpl + D {22) p p p» 2 + D^p^p" 

+ D§*> {p\p\ + pM + Df ] (p^ 3 + p»p%) + D$«> (f 2 pl + plpl) , 

DT vp = D^ig^pP + perm) + Df 2 \g p '"p p + perm) + Df 3 \g^p p 3 + perm) 

+ D { 3 ll2 \p p p v l p 2 + perm) + Dg 113 ^ (piPiP^ + perm) 
+ ^ 3 221) {p p 2 p 2 p{ + P erm ) + ^3 223) {P2P2P3 + P erm ) 

+ ^^(PsP&Pi + P erm ) + D 3 32) (P3P3P P 2 + perm) + D { 3 123) (p^p p 3 + perm) , 

(B5) 

where "+perm" indicates that the sum over all possible permutations of the tensor indices 
is understood. In the following we will give the full structure of the scalar box integrals, 
including both pole and finite parts, while for the corresponding tensor integrals we will 
only give the IR pole parts, since they can be of interest in checking the IR structure of the 
virtual cross section. We will write the pole part of each tensor integral coefficient as 

MA IR {Df l) ) , (B6) 

where J\f t is defined in Eq. ()17j) . and give for each box integral the non zero Am(D^), 
A IR (Di Jk) ), and A IR (D^> kl) ) coefficients. 



Df 


\JR- 


-pole 


I 

16n 2 




\lR- 


-pole 


i 

16vr 2 


D (JM) 


\IR- 


-pole 


i 

167T 2 



a. Box scalar integral DO (1,2) 

2,s 



The scalar integral appearing in diagram B^ s can be parameterized according to Eq. ()B3 
with: 



N 1 = k 2 , N 2 = (k + ptf-m 2 



t 1 



N 3 = (k + p t +p h ) 2 -m 2 t , N A = (k-p' t ) 2 -m 2 . (B7) 
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DO (2) is obtained from DO (v> by exchanging p t <-> p+. This integral also appears in the 
qq — ► tt/i calculation and has already been presented in Ref. 2JJ] . We repeat it here for 
completeness. 

The part of D0 R (i) which contributes to the amplitude squared is of the form: 

DO b^ = TTl^—f " -T7T, (— + Xo) , (B8) 

B 2, s I67H uj\ (cr — u)\ — 0J2 + Mfi) \ e ) 

where Af t is defined in Eq. (|17|). The pole part X_i is: 



where (3 t i is given in Eq. (J2EJ). The finite part X can be calculated using Ref. |45j. All 



tensor box integrals associated to i^ 1 ] and B^ s are IR finite. 



6. Box scalar integral DO (1,2) 

^7,4 



The scalar integral appearing in diagram B^ 1 ), D0 R (V), can be parameterized according 

' 7,t 



to Eq. (|B3I) with: 

N 1 = k 2 , N 2 = (k + qi f , 
N 3 = (k + gi - p 4 ) 2 - m 2 , iV 4 = (k + gi - p t - p fe ) 2 - m 2 . (BIO) 

The part of D0 R (V) which contributes to the virtual amplitude squared is of the form: 



7,4 



where the coefficients X_ 2 , A_i, and X are given by: 

al, = h. M , 
V win / 

X„ = i?e < -^ 2 + In 2 f ^) + In 2 - In 2 

2 1 u I V f^l + 2 In Vn ^ 



n / v^i/ v wi j v r i 

2 Li 2 ( Vl ~^~"* ) ~ 2 Li 2 ( " 1+ ; 2 ~ Tl ) + 2 Li 2 ( ^^-^ ) - Jo 



(B12) 



51 



with 



Io = In 



- Li 2 



mf 



+ 



A+(ti - r 2 ) + r 2 



-Li 2 
+ Li 2 



^2 

A+(ti - r 2 ) +r 2 



-r 2 - A + (ri - r 2 ) 

T\ - T 2 

+ (a+ - a_; 



(B13) 



and 



1 / / 4m? 



(B14) 



h 



The tensor integrals associated with B^l also contain IR divergences. Using the notation 



introduced in Eqs. ()B5|) and ()B6|) . only the following coefficients of Dl 



(i) 



f f f f 

+ 



2 ncJi e 2 ri^i 



-A T1 + 



T"2 



r 2 + Wi 



(Ar 2 -A W1 ; 



(B15) 



and of D2"V 



(ii), 

2 / 



111 1 

+ 



2 Tito>i e 2 tiuji 



A 



7 (t 2 + cO 2 



(A T2 -A W1 ) - 



T 2 + CJi 



(B16) 



are IR divergent. 



D0 B (2) and the corresponding tensor integrals are obtained from -D0 R (i) by exchanging 



U- 



qi <-> g 2 and p 4 <-> i.e. by exchanging n <-» r 2 and t^i <-> w 2 in Eqs. flBll|) - ()B16|) . 



c. Box scalar integral DO 



B. 



(1,2) 



The scalar box integral appearing in diagram B9~) , DO (i) , can be parameterized according 

' ''8,( 

to Eq. JHl with: 



iV 1 = A; 2 , N 2 = (k + gi ) 



N 3 = (k + qi + q 2 y , N 4 = (k + qi + q 2 - p' t ) 



l\2 



(B17) 



The part of D0 R (i) which contributes to the virtual amplitude squared is given by: 



D0 bW 



167P 



ar 2 



(B18) 
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where Mt is defined in Eq. (fT7f). and the coefficients X_ 2 , X_i, and X are given by: 
X_i = In 



3 
2 



a To 



Vm 2 



ln 2 (^ 



m 



2Li 2 ( 1 + — ) + — 

T 2 



7T 



(B19) 



The tensor integrals associated with B^] also contain IR divergences. Using the notation 



introduced in Eqs. (jB5|) and (jB6|) . only the following tensor coefficients of Dl (1 



B 



A^pf^ 



3 11 1 

2 CTT 2 e 2 (XT 2 
111 1 



K + A T2 + 



T 2 + UJi 



(A r2 - A 



2 (TT2 e 2 o~t 2 



A, 



(B20) 



of L>2*V 



(ll)v 
2 > 



a /r (^ 12) ; 



A /R ( J D 



(22) > 



3 1 1 



1 



T 2 



2 crr 2 e 2 ctt 2 [ r 2 + Co>i 
111 1 . 1 

2 + A -2- > 

o 2 + — (-1+A T2 )- 

2 crr 2 crr 2 e 



+ A a + A T2 + 



(r 2 + cui)' 



(A T2 -A^J 



(B21) 



and of L>3^ 



A/Hp. 



a„(^ 112) : 



(222) j 



3 11 1 

+ 



2 (XT 2 e 2 2ar 2 
2a; 2 



3r 2 



+ 



2r 2 o;i 



2A„ - 2A T 



r 2 + CJi 3 



72 ^'l [T 2 +U lj 
1 

e 



1 1 1 



J_ 1 

2 cxr 2 e z ar 2 e 

111 1 /1 A x 1 
2 + (!- A r 2 )- , 

2 <tt 2 e z ar 2 e 
111 1/3 
2^~? + ^ 2 \2~ T2 



(B22) 



are IR divergent. 

DO (2) as well as the corresponding tensor integrals can be obtained from D0 R (i) by 
exchanging q% «-> g 2 and p£ <-> i.e. by exchanging r x «-> r 2 and uj\ <-> cj 2 in Eqs. ()B18|) - 
flB22j) . 
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d. Box scalar integral DO n (i,2) 

The scalar box integral appearing in diagram Byj t , D0^ iot , can be parameterized accord- 
ing to Eq. IB3I with: 

N 1 = k 2 , N 2 = (k + gi ) 2 , 
N 3 = (k + q 1 -p / t ) 2 -m 2 , N A = (k + qi + q 2 -p' t ) 2 -m 2 . (B23) 

The part of DO w which contributes to the virtual amplitude squared is given by: 

^L^(^)(^ + ^ +A °) • (B24) 

where the coefficients X_ 2 , X-±, and Xq are given by: 

X-2 = \ , 



2 ' \ ™2 ] Xii \ ^2 



+ 21n(^±^)lnf ^ )+2ln T -^)ln f * 



T4 J \ r 2 + r 4 + ^l/ V T 2 J V T 2 + T 4 + ^1 

2 Li 3 [ V2 + T4 + ^ - 2 Li, f T2 + T4 + ^ - 2 Li ^ + ^ + ^ 



r 4 y \ r 2 / \ r 2 r 4 



(B25) 



The tensor integrals associated with 73^ t also contain IR divergences. Using the notation 



introduced in Eqs. (jB5|) and ()B6|) . the only IR divergent tensor coefficients of Dl (1 



o(l) • 
^10,4 



A /fi ( J DS 1) ) = ^ + -^—-{(r 2 + uj 1 )A Ti +u 1 (A T2 -A UJl )]- 

2 r 2 r 4 e z t 2 t 4 (t 2 + cji) e 



(B26) 



of D2» v m : 



10, t 

A IR (D ( 2 U) ) = i-Li + 1 [r 2 (r 2 + a*) - (r 2 + ^) 2 A T4 - (A T2 - A^)] - , 

(B27) 

and of D3 M ^ } : 

B 10,t 

A/*(^ 1U) ) = 4—4-5 ? T ^ [-2^(A T2 -A^)-2(r 2 + ^) 3 A, 4 

+ 3r 2 (r 2 + c^) 2 + 2r 2 w 1 (r 2 + wi)] - , (B28) 
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are IR divergent. 

DO (2) can be obtained from DO w by exchanging p' t <-> p t , i.e. by exchanging n <->• r 4 , 

"10,t "10,* 

r 2 <-> r 3 , and <-> tu 2 in Eqs (JB24J) and ()B28|) . 



2. Pentagon integrals 

The scalar and tensor pentagon integrals originating from the generic pentagon diagram 
Pi j in Fig. El are of the form: 

/' /! 1; ! , (B29) 



where 



jV 1 = (A; 2 -m 2 ) , iV 2 = (A; + pi) 2 -m 



1 5 



iV"3 = (fc + pi +p 2 ) 2 -ml , N A = (k + pi +p 2 +P3) 2 ~ ml , 
N 5 = {k + p 1 +p 2 +p 3 +p±) 2 -ml , (B30) 

Pi, P2, P3, Pa, and ps = —p\ — p 2 — Pz — P4 are the external (incoming) momenta connected 
to the pentagon topology, while m , mi, m 2 , m 3 , and m 4 are the masses of the propagators 
in the pentagon loop. 

The scalar pentagon integrals are evaluated as a linear combination of five scalar box 
integrals, using the technique originally proposed in Ref. 0,10. In particular, we use: 

1 5 

k=l 

(k) 

where each scalar box integral D0 P . can be obtained from the scalar pentagon integral 



E0p i . in Eq. (jB29J) by dropping one of the internal propagators. The coefficients c& are 
given by: 

5 

c k = J2 S h , (B32) 
1=1 

where Ski is the symmetric matrix: 

S kl ^ l -(Ml + Mf-p 2 kl ) , (B33) 

built out of the internal propagator masses M k and Mi and the linear combination of external 
momenta = pt + . . . + pf^ (k, I — 1, . . . , 5). A thorough explanation of this method is 



given in Ref. (21 



22, 



to which we refer for more details. 
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We write the tensor pentagon integrals as a linear combination of the linearly independent 
tensor structures built of the external momenta p±, p 2 , p%, and p%, which in d — 4 constitute 
a complete basis. Our notation for the pentagon tensor integrals is as follows: 

= #J V + E^pl + E {3) p» + E^pl , 
E7T = E^plpl + Ef 2) pW 2 + E^VsPs + ^pfrS 

+ E { 2 12) (plpl + p\p») + Ei 13) (jftp» + p*p£) + E^ ] [f xV \ + p\pl) 

+ e {23) (pgpg + pM) + 4 2A) (p%pZ + pU) + 4 U) WL + vlvl) , 

E3T> = E^ViPxP'x + Ef 22) p^ 2 p p 2 + EtMp^ + E^p'M 

+ E^ 12 ^ (P1P1P2 + perm) + E^ 113 ^ (piPiP^ + perm) + -E^ 114 '* (piPiPl + perm) 
+ E^ 21 ^ (p 2 p 2 Pi + perm) + E^ 23 ^ (p 2 p 2 p^ + perm) + Ef 24 ^ (p 2 p 2 p1 + perm) 
+ E^ 31 ^ (P3P3P1 + perm) + E^ 32 ^ (p^p^p 2 + perm) + Eg 334 -* (p^p^p^ + perm) 
+ Eg 441) + perm) + £^ 442) {p1PaP2 + perm) + E| 443) (p^plpl + perm) 

+ E^ 123 ^ (piP 2 Pz + perm) + E^ 2 ^ {piP 2 p1 + perm) + E^ 3 ^ (piP^pl + perm) 
+ E^ 234) (p^pIpI + perm) . 

(B34) 

The calculation of gg — > tih involves the six pentagon structures illustrated in Fig. |S| For 
each of them we will give in the following the IR pole parts of the corresponding scalar 

integrals, as well as the coefficient (in terms of the Ski matrix) and the IR singular box 

(k) 

scalar integrals D0 P . . out of which they can be calculated. We will moreover list the IR 
pole parts of the corresponding tensor integral coefficients, since they may be of interest in 
checking the IR structure of the virtual cross section. We will write the pole part of each 
tensor integral coefficient as 



Ehm-poie = ^ 2 NAm{E?) 

E?% R - po ie = -ir-MAmiE^) 



E^% R - pole = ^- 2 MAi R {Ef l) ) , (B35) 

where J\f t is defined in Eq. (II 7|) . and give for each pentagon integral the non zero Ajr(E^), 
Ai R (EV k) ), and A IR (E¥ M) ) coefficients. 

As in Section fB II we express our results in terms of the kinematic invariants a,Ti,u>i of 
Eqs. (JH and (jBTJ), and f3 tl of Eq. (HJ). 
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a. Pentagon scalar integral EOp 1 



The pentagon scalar integral arising from diagram P\ tt coincides with EQp 1 of the qq — > tth 



calculation of Ref. 



21 1, and can be parameterized according to Eq. (jB29|) with: 



N 2 = (k + qi ) 



l\2 



N 4 = (k + q l + q 2 - p' t ) 



m 



t i 



N 3 = (k + qi + q 2 ) 2 , 

N 5 = (k + qi + q 2 - Pt ~ Phf 



m. 



(B36) 



We summarize here for completeness the results obtained in Ref. j2l|. The Ck (k — 1, . . . , 5) 
coefficients in Eq. (|B31|) are obtained, according to Eq. (jB32|) . as: 



Ck 



(B37) 



i=i 



where 



S{Pi,t 



I 








—a 







\ 













t 2 


T\ 






—a 











—<jJ 2 






-u x 


r 2 





2m t 2 


(2m? - Ml) 




\ 





n 


— <jJ 2 


(2m? - Ml) 


2m? 


/ 



(B38) 



The part of E0p lt that contributes to the virtual amplitude squared can be written as: 

A_2 , X_i 



E0* 



16/T : 



Xn 



(B39) 



where X_ 2 , X_i, and X are obtained using Eqs. ()B31j) . ()B37|) . ()B38|) and the results for 

(k) 

the -DOp integrals presented in the following. The expressions for X_ 2 and X_i have the 
following form: 

X- 2 



1 2 
+ — 



1 

2a \ ^\ti nr-2 
1 



GTYT 2 

1 



(-A CT + A^ + K 2 - An - A T2 ) + 
(A T1 - A T2 + A W1 ) . 



(A T2 - A n + K 2 ) + 



(B40) 
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The tensor integrals associated with Pi t contain IR divergences. Using the notation intro- 



duce in Eqs. (jB34|) and (jB35|) . only the following coefficients of El p 



1 



1 2 \ 1 1 
1 - + - 



2<jt 2 \to 2 T\) e 2 
1 



+ 



(A 



T\ ' A T2 A^jy 



T\T 2 

1 

e 



( A(j + A T1 + A T2 - A W1 - A W2 ) 



(2), 



U 2 T 2 

11 1 / A A AN 1 

+ (A n - A T2 - K 2 ) - , 



of E2 P \-. 



&m(Ei 12) : 
A IR (Et y 



2ot 2 uj 2 £ 2 &t 2 uj 2 



2<JT 2 [^UJ 2 

^1 

t\t 2 {t 2 + U)\ 
1 1 



2ar 2 uj 2 e 2 
1 1 



2\ 1 1 


1 


T\ J e z a 




-(A W1 -A T2 ) 


1 




e 


(An - 


A T2 



(A, 



A T 



A, 



(JT 2 UJ 2 
1 



Au; 2 ) > 



2<tt 2 u; 2 e 2 ffr 2 w 2 



A r , + 



7l 



(r 2 + 



(A* 



A, 



(B41) 



T 2 UJ 2 



(A T2 + A W2 - A 7 



(B42) 



and of E3 P up : 



AzK(^ in) ) 



A^ 112) ) 
A ra (if 21) ) 

A^(^ 222) ] 



1 



2ar 2 
+ 



1 



2 

71 



1 



7"! T 2 (T 2 + Wl 



T"lT 2 

(V - A T 



(A 



^'2 



A 



n 



A CT ) + (A T2 + A W2 - A 



r 2 u-> 2 



1 1 /A 

+ (A T1 - A 



n(r 2 + wi) 
1 



2ar 2 uj 2 e 2 <jt 2 uj 2 
1 1 1 



rj " A W2 



1 1 



1 



2ar 2 u 2 e 2 cr 



r 2 o; 2 



A, 



-A„ + 



71 

(r 2 + ux] 



e 

(A. 



At 



r 2 o; 2 (ri + lu 2 )'- 



(A W2 - A T1 ) - 



r 2 (ri + lu 2 ) 



1 

— j 
e 

(B43) 



are IR divergent. 

(k) 

We present in the following the IR singular box scalar integrals D0 Pi , which are used in 
Eq. (|B31j) to calculate E0p l t . D0 p ^ t is finite and we will not discuss it further. 
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Box scalar integral D0^ t 
DOpj can be parameterized according to Eq. ()B3j) with: 

N 1 = k 2 , N 2 = (k + q 2 f , 
N 3 = (k + q 2 - p' t ) 2 - m 2 t , N 4 = (k + q 2 - p' t - p h ) 2 - m 2 . (B44) 

and can be obtained from DO a) in Section lE> 1 bl bv exchanging qi <-> q 2 and p t <-> p' t , i.e. 
by exchanging Ti <-> r 2 , and co>i <-> u; 2 . 
Box scalar integral DOp^ 



can be parameterized according to Eq. ()B3|) with: 

N 1 = k 2 , N 2 = (k + qi ) 2 , 
A^ 3 = (k + q 1 - p t ) 2 - m 2 , N 4 = {k + qi -p t - p h ) 2 - m 2 . (B45) 

and is equal to DO^m in Section IB 1 bl 

B 7,t 

Box scalar integral D0^ t 

DOpj can be parameterized according to Eq. ()B3|) with: 

Ni = k 2 , N 2 = (k + q 2 ) 2 , 
N 3 = (k + qi + q 2 ) 2 , N 4 = (k + q x + q 2 - p t f - m 2 . (B46) 

and is equal to D0^ gt in Section IB 1 cl 
Box scalar integral DOp 1 1 



can be parameterized according to Eq. (jB3|) with: 

N 1 = k 2 , N 2 = (k + gi ) 2 , 
N 3 = (k + gi + q 2 ) 2 , N A = (k + Ql + q 2 - p' t ) 2 - m 2 , (B47) 

and coincides with DO a) in Section IB 1 cl 
b. Pentagon scalar integral E0p 2t 

The pentagon scalar integral arising from diagram P 2>t can be parameterized according 
to Eq. (HSU with: 



N l = k 2 , N 2 = (k-p' t ) 2 -m 2 , N 3 = (k-p' t + q 2 ) 2 2 



m 



t j 



N 4 = (k - p' t + qx + q 2 ) 2 - m 2 , N 5 = (k - p' t + q x + q 2 - p h f - m 2 . (B48) 



59 



The Ck (k = l, . . . , 5) coefficients of Eq. (|B31|) are obtained, according to Eq. (jB32|) . as: 



(B49) 



1=1 



where 



S(P2,t 



( o 



r 2 





2m? 



r 2 -UJ\ 

,2 o™2 



2m 2 



2m 2 2mf — a 
2m 2 t 

uj\ 2m 2 — a 2m 2 
\ a x 



2m 2 
2m 2 





a i 
a 2 



2m 2 



Ml 



a 2 2m 2 - Ml 



(B50) 



2m 2 ) 



and we have defined 

«2 



i\2 



2m t - (p t + p' t ) 



2m 2 t - (qi - phf = 2m 2 t + u x - n 



2m\ - a + uji+uj 2 - Ml 



Ml 



(B51) 



The part of E0p 2t that contributes to the virtual amplitude squared can be written as: 



DO 



p 2 „ 



16tt 2 



+ x 



(B52) 



where X_i and Xq are obtained using Eqs. (|B31j) . (|B49)l - (|B51jl . and the results for the DQ^ k l 



Pi, 



integrals presented in the following. The expression for X-\ has the following form: 

1 



(B53) 



t 2 u 1 (a - ui - uj 2 + Ml) Pu ~ V 1 - Ptt. 

All tensor pentagon integrals associated with P 2t are IR finite. 

We present in the following the IR singular box scalar integrals DOp^ which are used in 

)® ' DO®, and DO® 



Eq. dESH) to calculate E0p 2 t . DO® , DO® , and DO® are finite and we will not discuss 



them further. 



(3) 



Box scalar integral DOp 2 
DO® can be parameterized according to Eq. (JB3)) with: 



(k-p' t + qi +q 2 ) 2 ■ 



k 2 

2 



No 



(k-Pt) -m. 



m 



t 1 



Ni = {k+p t f 



rn 



t ' 



(B54) 



and is equal to DO fi (i) in Section IB 1 al 



Box scalar integral DOp^ t . 
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DOp can be parameterized according to Eq. (jB3|) with: 

N 1 = k 2 , N 2 = (k-p' t ) 2 -m 2 t , 
N 3 = (k-p' t + q 2 ) 2 -m 2 t , N 4 = (k + p t ) 2 -m 2 , (B55) 



and can be written as 

D < = 16^ M ^ + X °) ' (B56) 
where the pole part X_i is given by: 



r 2 (a - wi - cu 2 + M 2 ) ^ V 1 - At 
while the finite part Xq can be found from Eq. (2.9) of Ref. ^| with the identifications: 



= m\ = m 2 — > m 2 , (B58) 



c. Pentagon scalar integral E0p. it 

The pentagon scalar integral arising from diagram P^ jt can be parameterized according 
to Eq. (IB29J) with: 



N 1 = k 2 



N 2 



N 4 = (k-p' t + q 2 -p h ) 2 -m 2 t , iV 5 = (fc + p t ) 2 -m 2 



(B59) 



The Cfc [k = 1, . . . , 5) coefficients of Eq. (|B31J) are obtained, according to Eq. (jB32|) . as: 



where 



1=1 







2m 2 



T-2 

2m 2 
2m 2 



t 2 2mf 

T\ a 3 2m 2 — 
\ a x a 2 



,2 itl 



\ 



Ml 



2m\ - a 2 
2m 2 

2m 2 2m 2 ) 



2m 2 



(B60) 



(B61) 
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and we have denned 



a 3 = 2m 2 - (q 2 - Ph) 2 = 2m 2 - M 2 + uj 2 + T\ - r 2 



(B62) 



while ai and a 2 are given in Eq. (jB51|) . 

The part of E0p 3 1 that contributes to the virtual amplitude squared can be written as 



EOp, 



P 3 ,, 



IQti 



+ X 



(B63) 



(fc) 



where X_i and X are obtained using Eqs. (jB31|) . ()B60|) - (|B62|) . and the results for D0 Pa 
given in the following. The expression for X_i has the following form: 

1 



1 , ( 1 + Ptt 
m 



nr 2 (a -u 1 -u 2 + M 2 ) (3 t i V 1 - Pu 



(B64) 



All tensor pentagon integrals associated with P 3t are IR finite. 

We present in the following the box scalar integrals D0^ t , which are used in Eq. (|B31j) 
to calculate E0p 3 t . DoQ t , D0^ t , and D0 [ ^ t are finite and we will not discuss them further. 

Box scalar integral DOp 3 1 



DOpzt can be parameterized according to Eq. ()B3|) with: 



Ni = k 2 



N 2 = (k — p t ) — m 



t j 



N 3 = (k - p' t + q 2 - p h ) 2 - m 2 , N A = (k + p t f - m 



2 

t 5 



(B65) 



and can be written as 



where the pole part X_ x is: 



16n 



1 



l. ln fl + fa 



(B66) 



n (a - ujx - u 2 + M 2 ) Pu V 1 - Ptt 
while the finite part Xq can be found from Eq. (2.9) of Ref. 4^ with the identifications: 



(B67) 



m = m 



m 



t j 



s ->• (Pt+ Pt) =a + m h -uj 1 -uj 2 
t -»• (qi ~ Pt) 2 = m 2 - n . 



(B68) 



Box scalar integral D0^ t . 
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_DOp 3 1 can be parameterized according to Eq. (jB3|) with: 

N 1 = k 2 , N 2 = (k-p' t ) 2 -m 2 t , 
N 3 = (k-p' t + q 2 ) 2 -m 2 t , N 4 = (k + p t ) 2 -m 2 , (B69) 

and is equal to -DOp in Section IB 2 bl 
d. Pentagon scalar integral E0p 4t 

The pentagon scalar and tensor integrals arising from diagram P^ t can be found from 
the corresponding integrals for diagram P 2 ,t by exchanging qy <-> q 2 and p t <-» pj, i.e. by 
exchanging r x <-> r 2 , r 3 <-> r 4 , and co>i «-> c<j 2 . 



e. Pentagon scalar integral E0p 5t 



The pentagon scalar integral arising from diagram P^ t can be parameterized according 



to Eq. (jEZSD with: 



Ni = k 2 , N 2 = (k + qi ) 



N 3 = (k + gi - p\Y - m 



« 5 



iV 4 = (fc + qi + <?2 - Pt) - m 2 , N 5 = (k + qi + q 2 -p' t - p h ) 2 - m 2 . (B70) 
The Cfc (k = l, . . . , 5) coefficients of Eq. (|B31|) are obtained, according to Eq. (jB32|) . as: 



1=1 



(B71) 



where 



( r 4 




r 2 





T\ 

a 3 



\ 



r 4 2m 2 2m 2 
-uj x t 2 2m 2 2m 2 t 2m 2 - Ml 
\ n a 3 2m 4 2 - Ml 2m 2 J 



(B72) 



with a 3 as defined in Eq. ()B62j) . 

The part of E0p 5t that contributes to the virtual amplitude squared can be written as: 



E0 



Ps,t 



167T 



+ l +X 



(B73) 
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where X_ 2 , X_i and X are obtained using Eqs. (jESH) , (IB7T1) . (IB721) . and the results for 
given below. The expressions for X_ 2 and X_i have the following form: 



2 



1 1 



2TiT 4 V UJi T 2 



— (A T2 - A n - A„J + — (A T2 + A T4 - A, 



(B74) 



TiT 4 

The tensor integrals associated with Ps ]t also contain IR divergences. Only the following 
tensor coefficients of El □ : 

Ajk(^ 1} ) = -^4 + ^^( A -i-A T2 -A r4 )i } 
£tit 2 t± e z T1T2T4 e 

(B75) 

ofE2- : 

A /K (4 11) ) = "^^^ + 7—— — r [( r 2 + ^i)A T4 — cji (A w — A r )] - , 



(B76) 



and of E3 H up : 



A - ( ^ m,) = 2^7- ~ nr 2 rl + ^ ^ + ^ + ^ + ^ ' B77 » 



(B78) 



are IR divergent. 

We present in the following the IR singular box scalar integrals D0 p ^ t , which are used 
in Eq. ()B31|) to calculate E0p 5t . D0p^ t and DOp*} t are finite and we will not discuss them 
further. 

(3) 

Box scalar integral D0 P ^ t 

DOp*} t can be parameterized according to Eq. (|B3|) with: 

iV\ = A; 2 , , N 2 = (k + gi ) 2 , 
N 3 = {k + qi + q 2 - p't) 2 - m 2 t , N A ={k + p t ) 2 - m 2 t , (B79) 

and coincides with DOp^ in Section fB 2 a[ after shifting k — > —k — q\. 
Box scalar integral DOp^ t . 
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D0^ t can be parameterized according to Eq. (jB3|) with: 



N\ = k , N 2 = (k + qi ) 2 , 
N 3 = (k + gi - p' t ) 2 - m 2 t , N A = (k + p t f - m 2 t . (B80) 
The part of -DOp which contributes to the virtual amplitude squared is given by: 

<!.-I^,(^)(^ + ^ + *„) , (BSD 

where the coefficients X_ 2 , X_ 1: and X are given by: 

X- 2 = 1 , 

X_ 1 = -ln^-lnf^ 

ml) \mf 



with 



/. Pentagon scalar integral E0p ( 



a.t 



(B82) 



z ± =-(l±A) , A^l-^f-, (B83) 
and 03 defined in Eq. flB62|) . 

(5) 

Box integral .DOp^. 

DOp^ can be parameterized according to Eq. (jB3|) with: 

iV 1 = fc 2 , N 2 = (k + gi ) 2 , 
N 3 = (k + qi -p' t ) 2 -m 2 t , N A = (k + qi + q2 -p' t ) 2 -m 2 t , (B84) 

and is equal to D0 R (i) in Section IB 1 dl 



The pentagon scalar integral arising from diagram P 6 t can be parameterized according 
to Eq. (|B"29j) with: 



N 1 = k 2 , N 2 = (k + qi y , N 3 = (k + qi -p' t y -m 



t 5 



N 4 = (k + qi - p' t - p h ) 2 - m\ , N 5 = (k + q 1 -p' t -p h + g 2 ) 2 - m 2 . (B85) 

We note that E0p 6t can be obtained from E0p 5t by shifting k — > —k — qi and exchanging 
p t <-» pj, or equivalently by exchanging ri <-» r 4 , r 2 <-> r 3 , and co>i <-> c<j 2 . The same applies to 
the tensor pentagon integrals El p , E2 p u 6t , and E3 p u 6 P t . 
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APPENDIX C: PHASE SPACE INTEGRALS FOR THE EMISSION OF A SOFT 
GLUON IN THE TWO CUTOFF PSS METHOD. 



In this Appendix we collect the phase space integrals for a final state soft gluon that 
are used in calculating the results reported in Eq. (|39|). We parameterize the soft gluon 
(i-momentum in the gg rest frame as: 

k = E g (l, . . . ,sin9 1 sm9 2 ,sm9 1 cos92,cos9i) , (CI) 

such that the phase space of the soft gluon in d = 4 — 2e dimensions can be written as: 



d{PSg)soft 



T(l - e) tt 6 
T(l - 2e) (2vr)3 J, 



<5 sx /i/2 



dE a E}~» x 



d0i sin 1-2e 0i / d9 2 sm~ 2e 9 2 . 



(C2) 



Then all the integrals we need are the following four: 



d(PS, 



'9) soft 



(qvk)(q 2 -k) (4tt) 



;M2 



- -ln(5 s ) - -A c 

e 



e 2 e 



d(PSg) SO f t 



(qi-pt) 



-^ + i(A 2 +4A CT ln(5 s )+41n 2 (<y) 
1 



(q v k)( P fk) ' (4tt) 



1 2 A 

-7 - -A, 



e e 3 



^A 2 + 2A T1 K + 2 ln 2 (5 s ) + 4A ri ln(5 s ) + F(q llPt ) 



d(PSg) SO ft 



(PfP't) 



(prkM-k) ' ~ (4tt) 



s« - 2m 2 



— + 2A (7 + 41n(5 s 
e 



■A tf - 



4 



A 2 f - — Li 2 

Ptt U Ptt ' 



2A 



1 + Ptt 



d(PS 6 



Pt 



1 



' 9>sof \ Pt -kf " ' (4tt) 2 
where s t t, Ptt and A t j are defined in Eq 
function: 



-- + 2A CT + 41n(<5 s )-2-j-A tt - 
e Pit 



(C3) 



= In 2 ( : 

\1 — Pf cos9if 



Moreover we have denoted by F(pi,pf) the 

i + /V 



-In 



+2Li 2 



^(l-cos^; 



-2Li 2 



/?/(! + cos%; 



(C4) 



1 — Pf J \ 1 — Pf cos Oif 

where cos0j/ is the angle between partons i and / in the center-of-mass frame of the initial 

state partons, and 

m? „ Sif 



Pf 



(P° f ) 2 



1 — Pf cos 9i f 



P° f Vs 



(C5) 
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All the quantities in Eq. (|C4|) can be expressed in terms of kinematical invariants, once we 
use Sif = 2pi-pf and: 



Pt 



o s-sth + m; , o s-s t h + m% 
and Pi = 



2^ 



2^ 



(C6) 



with s fh = (pf +p h ) 2 - 



APPENDIX D: INTEGRATED SOFT FUNCTIONS FOR THE ONE-CUTOFF 
PHASE SPACE SLICING METHOD. 

In this appendix we give the explicit form of the integrated soft functions S a b in the 
three possible cases in which: both partons (a, b) are massless, one is massless and the other 
is massive, and when both are massive. These expressions have been originally presented 



in Refs. |2 
h — > qqti + g in Ref . 



28 1, and used in the calculation of the soft part of the real cross section for 



2l|. 



When both partons a and b are massless S a b is simply given by 



Sat = -T-N- 



Sab 



1 



(Dl) 



2tT £j y S m i n J \S m in J ^ 

For h — > ggti + (g, q, q), this occurs when a and b correspond to the two hard gluons g A and 
g B of Eq. (|54|) . in which case s a t, = s is the partonic center of mass energy. 

When one parton is massive and the other is massless, the function Sab has the form : 



S ab = 



2tt rn 



47T/i 2 



1 / s ab 



mi 



Sab 
, Smin 
1 f Sgb 

2e \m 2 t 



x 



m. 



Ots « T I 

—N 

2tt T(l - e) 



47T/i 2 



1 1 

2^2 + 2~e 



Ilnf^ 



+ Tin 2 



Safe 



-In" + t In + - In 



Safe 



Safe 



:d2) 



For /i — > ggtt+ (g, q, q), this occurs when a— 1,2 (where 1, 2 denote the initial gluons g A , g B ) 
and b — t,t, and we therefore have four possible integrated soft functions of this type: S u , 
Sit, S 2t , and S 2 t- 
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Finally, when both partons are massive, i.e. when a = t and b = t, the integrated soft 
function S t t is 1 281: 



27T m-e I S min I J \++ \ 6 



where we have defined: 



--J, = 1 - 



Stt 



J a = J s ln 



(2m 2 t + sdfr 



r 2 + Xu 

2 I ' 

Smin^t 



-A tt - 



J b = (r+ - r_) [1 - 2 ln(r+ - r_) - ln(r 4 



+ 



Li 2 ( 1--^ I -Li 2 l 1-- 



+ (r- + r+ 

Ptt and A 4 f are defined in Eq. (|2~Hj) while A# and r± are given by: 



+ 1 + T_T_ 
1 



l-ln(r + )ln(r_) + ^ln 2 (r_ 



(D3) 



(D4) 



A tt - 



4m 



2m? 



± 



2m? 



(D5) 
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